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1 Introduction 

For the plasma physicist, the solar corona offers an outstanding example of 
a space plasma, and surely one that deserves a lifetime of study. Not only 
can we observe the solar corona on scales of a few hundred kilometers and 
monitor its changes in the course of seconds to minutes, we also have a wide 
range of detailed diagnostics at our disposal that provide immediate access 
to the prevalent physical processes. 

Yet, solar physics offers a rich field of unsolved plasma-physical problems. 
How is the coronal plasma continuously heated to > 10 6 K? How and where 
are high-energy particles episodically accelerated? What is the internal dy- 
namics of plasma in magnetic loops? What is the initial trigger of a coronal 
flare? How does the corona link to the solar wind, and how and where is the 
latter accelerated? How is plasma transported into the solar corona? 

Why, then, study stellar coronae that remain spatially unresolved in X- 
rays and are only marginally resolved at radio wavelengths, objects that 
require exposure times of several hours before approximate measurements 
of the ensemble of plasma structures can be obtained? 

There are many reasons. In the context of the solar-stellar connection, 
stellar X-ray astronomy has introduced a range of stellar rotation periods, 
gravities, masses, and ages into the debate on the magnetic dynamo. Coronal 
magnetic structures and heating mechanisms may vary together with varia- 
tions of these parameters. Parameter studies could provide valuable insight 
for constraining relevant theory. Different topologies and sizes of magnetic 
field structures lead to different wind mass-loss rates, and this will regulate 
the stellar spin-down rates differently. 

Including stars into the big picture of coronal research has also widened 
our view of coronal plasma physics. While solar coronal plasma resides typi- 
cally at (1 — 5) x 10 6 K with temporary excursions to « 20 MK during large 
flares, much higher temperatures are found on some active stars, with steady 
plasma temperatures of several tens of MK and flare peaks beyond 100 MK. 
Energy release in stellar flares involves up to 10 5 times more thermal en- 
ergy than in solar flares, and pressures that are not encountered in the solar 
corona. 
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Table 1. Symbols and units used throughout the text 



Symbol, acronym Explanation 



R t or R 


Stellar radius [cm] 


Rq 


Solar radius [7 x 10 1() cm] 


c 


Speed of light [3 x 10 10 cm s^ 1 ] 


k 


Boltzmann constant [1.38 x 1CT 16 erg cm 3 K _1 ] 




Mass of electron [9.1 x 1CT 28 g] 


d 


Distance [pc] 


P 


Pressure [dyne cm" 2 ] 


L 


Coronal loop semi-length [cm] 


T 


Coronal electron temperature [K] 


T b 


Radio brightness temperature [K] 


n,N 


non-thermal electron density (N: integrated in energy) 


N 


Number of flares 


n e 


Electron density [cm -3 ] 




Hydrogen density [cm -3 ] 


B 


Magnetic field strength [G] 


f 


Surface filling factor [%] 


r 


Loop area expansion factor (apex to base) 


E 


Energy [erg] 


L R 


Radio luminosity [erg s -1 Hz -1 ] 


Lx 


X-ray luminosity [erg s -1 ] 


Lbo\ 


Stellar bolometric luminosity [erg s -1 ] 


A — Aol A 


Cooling function [erg s -1 cm 3 ] 


r 


(Decay) time scales, also: optical depth 


V 


Radio frequency 


Vp, UJp 


Plasma frequency, angular plasma frequency 


Vc, He 


Gyrofrequency, angular gyrofrequency 


K 


Absorption coefficient 


7 


Lorentz factor 


L B 


Magnetic scale height 


HRD 


Hertzsprung-Russell Diagram 


EM 


Emission Measure 


Q, DEM 


Differential Emission Measure Distribution 


EMD 


(discretized, binned) Emission Measure Distribution 



The present chapter provides a "stellar astronomer's view" of magnetic 
coronae. It cannot replace the knowledge of detailed physical processes in the 
solar corona, nor do we expect to find entirely new physical concepts from 
the study of stellar coronae without guidance by solar physics. However, as I 
hope to show in the following sections, stellar astronomy has provided some 
unexpected and systematic trends that may well help understand systematics 
of coronal behavior across a large range of stellar parameters. In this sense, 
the goal of this chapter is to present an overview of the basic, observed stellar- 
coronal phenomena and their interpretation analogous to solar physics, rather 



Stellar Coronae 3 



than the derivation of the basic plasma-physical mechanisms themselves. For 
a deeper understanding of the latter, I must refer the reader to appropriate 
textbooks and lectures on solar coronal physics. 

2 Stellar coronae - Defining the theme 

We consider coronae to be the ensemble of closed magnetic structures above 
the stellar photosphere and chromosphere together with their plasma content, 
regardless of whether the latter is the thermal bulk plasma or a non-thermal 
population of accelerated particles. Although often narrowed down to some 
specific energy ranges, coronal emission is intrinsically a multi-wavelength 
phenomenon revealing itself from the mctcr-wave radio range to gamma rays. 
The most important wavelength regions from which we have learned diagnos- 
tically on stellar coronae include the radio (decimetric to centimetric) range 
and the X-ray domain. The former is sensitive to accelerated electrons in mag- 
netic fields, and that has provided the only direct means of imaging stellar 
coronal structure, through very long baseline interferometry. 

The soft X-ray (0.1-10 keV) diagnostic power has been instrumental for 
our understanding of physical processes in the hot, magnetically trapped 
coronal plasma, and the recent advent of high-resolution X-ray spectroscopy 
with the Chandra and XMM-Newton X-ray observatories is now accessing 
physical parameters of coronal plasma directly. The adjacent extreme ultra- 
violet (EUV) range contains diagnostic relevant for the same temperature 
range as X-rays. 

3 The coronal Hertzsprung-Russell diagram 

Before discussing specific physical problems in stellar coronal physics, I will 
briefly review the phenomenology of stellar radio (here: 1-20 GHz) and X- 
ray emission and summarize stellar classes that are prolific coronal emitters. 
Fig-fflpresents Hertzsprung-Russell diagrams (HRD) of detected X-ray (left) 
and radio (right) stars. They show all basic features that we know from an 
optical HRD. Although the samples used for these figures are in no way 
"complete" (in volume or brightness), the main sequence is clearly evident, 
and so is the giant branch. The cool half of the subgiant and giant area is 
dominated by the X-ray and radio-strong sample of RS CVn and Algol-type 
close binaries. The top right part of the diagram, comprising cool giants, is 
almost devoid of X-ray detections, however (although well populated by radio 
em itters). The so-called c orona vs. wind dividing line (dashed in Fig. QJi; af- 
ter lLinskv fc Haischll97d) separates coronal giants and supergiants to its left 
from stars with massive winds to its right. It is unknown whether the wind 
giants possess magnetically structured coronae at the base of their winds 
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- the X-rays may simply be absorbed by the overlying wind material. Ad- 
ditionally, a very prominent population of (presumably) coronal radio and 
X-ray sources just above the main sequence is made up of various classes of 
pre-main sequence stars, such as classical and weak-lined T Tauri stars. 

Moving toward A-type stars on the HRD, one expects, and finds, a sig- 
nificant drop of coronal emission owing to the absence of magnetic dynamo 
action in these stars. However, this is also the region of the chemically pecu- 
liar Ap/Bp stars that possess strong magnetic fields and many of which are 
now known to be non-thermal radio sources as well. Finally, the very lumi- 
nous radio and X-ray emissions from O and B stars are believed to originate 
in non-magnetic stellar winds - we will not further discuss these stars. 




Fig. 1. Left: Hertzsprung- Russell diagra m based on about 2000 X-ray detected 
stars extracted from survey catalogs (see|Gudcl 2004] for references). The size of 
the circles characterizes logLx as indicated in the panel at lower left. The ranges 
for the spectral classes are given at the top (upper row for supergiants, lower row 
for giants), and at the bottom of the figure (for main-sequence stars). - R ight: 
Similar, but for 440 radio stars detected between 1-10 GHz (after iGiidefeooai 



4 Nonflaring radio emission from stellar coronae 

"Quiescent" (non-flaring) radio emission at levels of 10 12 — 10 16 erg s _1 Hz -1 
from magnetically active stars was entirely unanticipated but constitutes an 
important achievement in stellar radio astronomy: There simply is no solar 
counterpart! Quiescent emission can be defined by the absence of impulsive, 
rapidly variable flare-like events. Common characteristics of quiescent emis- 
sion are (i) slow variations on time scales of hours and days, (ii) broad-band 
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radio spectra, (iii) brightness temperatures in excess of coronal temperatures 
measured in X-rays, and usually (iv) low polarization degrees. 

4.1 Bremsstrahlung 

The Sun emits steady, full-disk, optically thick thermal radio emission at 
chromospheric and transition region levels of a few 10 4 K. However, such 
emission cannot be detected with present-day facilities, except for radiation 
from the very nearest stars, or giants subtending a large solid angle. Using 
the Rayleigh- Jeans approximation for the flux density S 

_ 2kTTv 2 irR 2 ( T \/ v \* f R\ 2 flpcV 

(i) 

(R = source radius, T = electron temperature, d = stellar distance, k = 
Boltzmann constant, v = observing frequency, r = optical depth), we find 
for optically thick chromospheric emission (with r = l,T = 1.5 x 10 4 K, 
v = 8.4 GHz) 

„ 0.05 / R \ 2 T 

Optically thin free-free emission from the hot, X-ray emitting plasma in 
coronal loops can be estimated as follows: The radio optical depth is 

idl ~ / n e dl ( 3 ) 

and the X-ray volume emission measure EM, using a filling factor / in the 
approximation of a small coronal height, 

EM X = 2R 2 irf J n 2 e dl. (4) 

(the factor of two accounts for the visible and invisible hemispheres, assuming 
that the EM is uniformly distributed). For plasma between 1-20 MK, the EM 
can be estimated fr om the X-ray luminosity Lx with a conversion factor A(T) 
ijMewe et al.lll985l A larger for lower T) 

L x = A(T)EM X re (1.5 - 6) x 10" 23 EM X [erg s" 1 ] (5) 

Thus, we find (using T MK = T/[10 6 K]) 

S = 2.6 x 10~ 52 ^ mJy ps (4-17) x 10- 30 ^^ mJy. (6) 

A(T)T MK fd 2 c ^mk/^pc 



Coronal bremsstrahlung contributions are presently out of reach for almost 
all stars. 
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4.2 Gyroresonance emission 

Because active stars show high coronal temperatures and large magnetic fill- 
ing factors that prevent magnetic fields from strongly diverging with increas- 
ing height, the radio optical depth can become significant at coronal levels 
owing to gyroresonance absorption. This type of emission is ob served above 
solar sunspots. The optical depth for the sth harmonic is (after IWhite et all 
1994) 

T[S,V) 4 v 2 2 *- 2 s![(s- l)!] 2 \2fi 
= 1.45^^1-^(1.3 x 1(T 4 ,4.1 x HT 6 .1.8 x 1(T 7 ) ( — ^— Y ( 7 ) 

where the three coefficients in the parenthesis are for the harmonics s — 3, 4, 
and 5, respectively, of the gyrofrequency, ^ghz is the observing frequency in 
GHz, n e is the electron density of the emitting hot plasma at temperature 
T, and \x = m e c 2 /kT. The magnetic scale height, Lb, is not precisely known 
but will be assumed to be 0.3i?* (see arguments in I White et alJllQM) . Using 
typical X-ray derived temperatures ([1—3] x 10 7 ) and X-ray emission measures 
EM = n 2 V (to estimate electron densities in the emitting volume), one finds 
that t invariably reaches unity at s — 3, 4, or 5 for an observing frequency of 
15 GHz while t < 1 for larger s. The highest harmonic that is still optically 
thick is relevant. 

The measured flux from the optically thick layer is then again given by 
CJ. Observations have shown the following: 

The observed fluxes at 15 GHz from M dwarfs are not compatible with 
the above prediction if the emitting layer of hot plasma covers the en- 
tire stellar surface. The observed flux is typically much smaller. The mi- 
crowave spectra ar e falling, not compatible with optically thick radiation 
l|White et alJll994l - see Fig.0). 

There are exceptions in which a rising spectrum from 5 to 15 GHz (Fig.[SJ) 
could be explaine d by gyroresonance em ission, while the lower-frequency 
spectrum cannot l|Giidel fc Benzlll989h . The same stars may show this 
feature only temporarily. 

These results imply that the hot plasma in general is not coincident with the 
strong-magnetic field regions in the corona. It must reside in lower- B regions 
in between or above the cooler regions, possibly implying rather extended 
coronal structures seen in X-rays. This plasma could be induced by flares in 
which magnetic loops reconnect to larger structures in w hich the low-densit y 
gas will not produce appreciable gyroresonance emission jWhite et alJll994j) . 
The cooler plasma also usually observed in active stars might be trapped in 
the strong fields but its gyroresonance emission is also negligible. 
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Frequency (GHz) 

Fig. 2. Radio spectra of the RS CVn binary HR f099 (upper set) and of the 
dMe dwarf UV Cet (lower set) at different flux levels. The gently bent spectra are 
indicative of gyrosynchrotron emission, and the high-frequency part of U-shaped 
spectra for UV Cet has been interpreted as a gyroresonance component (HR 1099 
spectra: courtesy of S. M. White). 



Lower-frequency radio emission cannot be due to gyroresonance absorp- 
tion: The radius of the optically thick layer would still at s = 3 — 5, but 
it would then be located at > 3i?„ (for dMe stars). If we extrapolate the 
the corresponding magnetic fields of more than 100 G down to photosphcric 
levels, we would find photosph eric field strengths much in excess of those 
observed ijGarv fc Lmskvlll98llh 

4.3 Gyrosynchrotron emission 

We could, however, allow for much higher T c g. The optically thick layer 
would in that case shift to harmonics above 10, the range of gyrosynchrotron 
radiation. The optica lly thick source sizes are then more reasonable for M 
dwarfs, with R i?* ( Linskv fc Garvlll983l) . and the optically thin emission 



may still be strong enough for detection. However, for a thermal plasma, the 
spectral power drops like v~ s at high frequencies, in contrast to observed 
microwave spectra that show ^-(P- 3 --- 1 ) for magnetically active stars (Fig. 01. 

Instead, the electron population could be non-thermal. Electron energy 
distributions in cosmic sources are often found to f ollow a power law, and 
this also holds for solar and stellar microwave flares l|Kundu fc Shevgaonkarl 
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n(j)=K( 1 -l)- s (8) 

There is wide support for this model from estimates of the brightness tem- 
perature (e.g., based on the stellar radius, or from resolved interferometric 
images) . The important question then is how these coronae are continuously 
replenished with high-energy electrons. 

The radio-spectral time development implied by this model can be analyt- 
ically calculated for the case of a sho rt injection of electrons into the corona 
ijChiuderi Drago fc Franciosinill 1993^1 . The outline of the derivation is as fol- 
lows: For trapped electrons with nonthermal number density n, an equation 
of continuity applies in energy, 



dn(j,t) d 
dt d"f 

with the solution 



(9) 



n(^,t)dj = n(7o)c£yo- (10) 

Here, n(7o)c?7o is the initial distribution. We now need an expression for the 
energy loss rate of an electron, d-f/dt. In an extended magnetosphere, the 
most relevant loss mechanisms are synchrotron loss and Coulomb-collisional 
loss, given by, respectively, 

- 7coii = 5 x l(T 13 n e [s- 1 ], 7-coii = 2 x 10 12 ^- [s] (11) 

n e 

-7B = 1.5xl0- 9 BV [b- 1 ], t b = 6 - 7 r X 2 1Q8 [s] (12) 

(where n e is the ambient thermal electron density: see, e.g.. IPetrosianl [l985: 
we have assumed a pitch angle of 7r/3). The total energy loss rate is thus 

^=a + /3 7 2 (13) 

with the appropriate coefficients from (|11|) and l|12fl . This equation has an 
analytical solution for n(-y,t) if the initial distribution is a power law as given 
by (|HJ), namely 



2/1 1^7(1 + At&nAat) - A + tan Aat}~ 

(1 - A-y ta,n Aat) 2 



n( 7 , t) = K(l + tan 2 Aat)A> ^ ^ Z 2. ^ l^T J ( 14 ) 



where A — ((3/a) 1 / 2 , and the initial power-law has been bounded by 70.1 < 70 
(typically a low energy, e.g., 70,1 = 1.1) and 70,2 > 7o, typically a very high 
energy, formally including 70.2 = 00. The evolution of the initial power- law 
boundaries can be computed; after a finite time, all non-thermal electrons 
have thermalized. Some characteristic results of these calculations are shown 
in Fig. El 
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Figure shows the flattening with time of the electron energy distri- 
bution due to the collisional losses at lower energies, and a cut-off due to 
synchrotron losses at higher energies. Figure [3J3 illustrates the intensity of 
radiation after calculation of the emissivity and the absorption coefficient, 
together with an observed spectrum of the RS CVn binary UX Ari. In the 
early spectrum, the optically thick portion is well visible on the low-frequency 
side. As the electron energy decays, the spectrum becomes optically thin, and 
once the high-energy electron population is depleted, the spectrum begins to 
fall off steeply as there are no electrons left emitting at the observing fre- 
quency. 




Fig. 3. Time-dependent spectral model of an RS CVn magnetosphere after injec- 
tion of a power-law electron population. Left: Evolution of the electron energy distri- 
bution due to collisional and synchrotron losses, starting from a power law. Middle: 
Evolution of the microwave spectrum with time, for a magnetic field of B = 10 G 
and an initial electron power-law index of 5 = 2. The electron density in the thermal 
power- law is n e = 2 x 10 8 cm~ 3 . Right: Time evoluti on of the radiation intensity at 
5 GH z, for different magnetic field strength. (From Ichiuderi Draeo fc Franciosinil 
Il993l . courtesy of E. Franciosini.) 



The time development of the emitted intensity as a function of mag- 
netic field strength (Fig. shows that the initial radiation originates in 
fields of higher strength (the "core"), but the responsible electron popula- 
tion decays rapidly due to synchrotron losses. Eventually, radiation from the 
slowly-decaying population in the weaker, more extended magnetic fields (the 
"halo" ) begins to dominate on longer time scales. 

A power-law index of 5 = 2 — 3 (as also used above) appears to usually 
fit active-stellar m i crowave spectra su ch as those shown in Fig. quite well 



ilJll998l) . This suggests ver y hard electron dis - 
tributions, similar to those seen in gradual solar flares l|Cliver et IDE986). 
But such distributions are observed during so-called quiescence - an indi- 
cation that quiescent radio emission is due to a flare-like process? We will 
return to this hypothesis in later sections. 
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5 Thermal X-ray emission from stellar coronae 
5.1 High-resolution X-ray spectroscopy 

The high-resolution X-ray spectrometers on XMM-Newton and Chandra 
cover a large range of spectral lines that are temperature- and partly density- 
sensitive. The spectra thus contain the features required for deriving X-ray 
emission measure distributions, abundances, coronal densities, and opacities. 




Fig. 4. Three high-resolution X-ray spectra of stars with largely differing activity 
levels: HR 1099, Capella, and Procyon. Data from XMM-Newton RGS. 



Figure 0] shows examples of X-ray spectra. The stars cover the entire 
range of stellar activity: HR 1099 representing a very active RS CVn system, 
Capella an intermediately active binary, and Procyon an inactive F dwarf. 
The spectrum of HR 1099 reveals a considerable amount of continuum and 
comparatively weak lines, which is a consequence of the very hot plasma in 
this corona (T w 5 — 30 MK). Note also the unusually strong Neix/Fexvn 
and Nex/Fexvn flux ratios if compared to the other stellar spectra (these 
ratios are due to an intrinsic compositional anomaly of the HR 1099 corona). 
The spectrum of Capella is dominated by Fexvn and Fexvin lines which are 
preferentially formed in this corona's plasma at T w 6 MK (Fig. [5J. Procyon, 
in contrast, shows essentially no continuum and only very weak lines of Fe. 
Its spectrum is dominated by the H- and He-like transitions of C, N, and 
O formed around 1-4 MK. The flux ratios between H- to He-like transitions 
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are also convenient temperature indicators: The Ovin A18.97/Ovn A21.6 
flux ratio, for example, is very large for HR 1099 but drops below unity for 
Procyon. 



Emissivities of Fe and O 

0.5 F 1 1 1 1 1 1 1 1 1 ' 1 1 ' 




log(Temperature) 

Fig. 5. Emissivities of several X-ray transitions, in particular of highly-ionized iron 
and oxygen species (Fig. courtesy of A. Telleschi, after Telleschi et al. 2005). 



5.2 Thermal coronal components 

The large range of temperatures measured in stellar coronae has been a chal- 
lenge for theoretical interpretation. Whereas much of the solar coronal plasma 
can be well described by a component of a few million degrees, magnetically 
active stars have consistently shown a wide distribution of electron temper- 
atures, reaching values as high as 50 MK outside obvious flares. 

The flux <f>j observed in a spectral line from a given atomic transition can 
be written as 



And 2 



(15) 



where d denotes the distance, and Gj(T) is the "line cooling function" (lumi- 
nosity per unit EM; Fig. that contains the atomic physics of the transition 
as well as the ionization fraction for the ionization stage in question, and A is 
the abundance of the element with respect to some basic tabulation used for 
Gj, For a fully ionized plasma with cosmic abundances, the hydrogen density 
Tin ~ 0.85n e . The expression 

Q(T) = (16) 

defines the differential emission measure distribution (DEM). I will use this 
definition throughout but note that some authors define Q'(T) = n e nndV/dT 
which is smaller by one power of T . For a plane-parallel atmosphere with 
surface area S, Eq. (fTTTf) implies 
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l 



Q(T) = n e n H SH{T), 



(17) 



where H is the temperature scale height. 
5.3 Observational results 

"Discrete" (binned) emission measure distributions (EMDs) reflecting the 
full DEMs can be obtained by various inversion algorithms from the observed 
spectra. Most EMDs have generally been found to be singly o r doubly peaked 
and c onfined on either side approximately by power laws (e.g.. lSchrii ver et alJ 
119951) . Interestingly EMDs are often very steep on the low-T side, and this 
is particularly true for the more active stars. For example, the slope of the 
stellar EMD in Fig. follows approximately Q cx T 3 on the low-T side. 

It is notable that the complete EMD shifts to higher temperatures with 
increasing stellar activity as seen in Fig. HJic, often leaving very little EM at 
modest temperatures and correspondingly weak spectral lines from ions of C, 
N, and O. 

As a consequence, a relatively tight correlation between the characteristic 
coronal temperature (e.g., derived from a DEM) and the normalized coronal 
luminosity Tx/Tbol is found: Stars at higher activity levels support hotter 
coronae, with the most active stellar coronae reaching characteristic temper- 
atures of several tens of MK. An example of solar analogs is shown together 
with the Sun itself during its activity maximum and minimum in Fig.0 Here, 



where Lx denotes the total X-ray luminosity (but EM and T refer to the 
"hotter" component in standard 2-T fits to ROSAT data). 

5.4 Interpretation of differential emission measure distributions 

Equations (|15|) and (|16|) introduce the DEM as the basic interface between 
the stellar X-ray observation and the model interpretation of the thermal 
source. It contains information on the plasma temperature and the density- 
weighted plasma mass that emits X-rays at any given temperature. Although 
a DEM is often a highly degenerate description of a complex real corona, 
it provides important constraints on heating theories and on the range of 
coronal structures that it may describe. Solar DEMs can, similarly to the 
stellar cases, often be approximated by two power laws Q(T) oc T s on either 
side of their peaks (Fig.[|)|). The Sun has indeed given considerable guidance 
in physically interpreting the observed stellar DEMs. 



Lx oc T 



i4.5±0.3 



(18) 
(19) 



EM cx T 



>5.4±0.6 
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log T 



Fig. 6. Left: Emission measure distributions of two intermediately active stars (bul- 
lets) and the Sun (histograms; EM is in units of 10 50 cm -3 ). The histograms refer 
to full-disk solar EMDs derived from Yohkoh images at solar maximum, including 
also two versions for different lowe r cutoffs for the in tensities in Yohkoh images, 
(figure courtesy of J. Drake, after iDrake et alJ 1200(1) . - Upper right: Calculated 
differential emission measures of individual static loops. The solid curves refer to 
uniform heating along the loop and some fixed footpoint heating flux, for differ- 
ent loop half-lengths labeled above the figure panel in megameters. Th e dashed 
curves illustrate the analytical solutions presented bv lRosner et alJ l)l97ST) for uni- 
form heating. The dotted lines show solutions ass uming a heating scale height o f 
2 x 10 9 cm (figure courtesy of K. Schrijver, after ISchriiver fc Aschwandenl l2002"). 
- Lower right: Examples of discrete stellar emission measure distributions of solar 
analogs (figure courtesy of A. Telleschi). 



The DEM of a static loop 

The DEM of the plasma contained in a static ma gnetic loop follows f rom the 
hydrostatic equilibrium (see, e.g., derivation bv lRosner et alJ|l978|) . Under 
the conditions of negligible gravity, i.e., constant pressure in the entire loop, 
and negligible thermal conduction at the footpoints, 

Q(T) oc p T 3 / 4 -*/ 2 +«- — (20) 

(l - [T/T a ] 2 -^) ' 



14 Manuel Giidel 





solar analogs 
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Fig. 7. Coron al temperatur e vs. X-ray luminosity for solar analogs. For details and 
references, see Giidel (2004). 



et alJll99lj) where T a is the loop apex temperature, and a and f3 are 
power-law indices of, respectively, the loop cross section area S and the heat- 
ing power q as a function of T: S(T) = SoT a , q(T) = qoT@, and % is the expo- 
nent in the cooling function over the relevant temperature range: A(T) cx T x . 
If T is not close to T a and the loops have constant cross section (a = 0), we 
have 

Q(T) cx T 3 / 4 -*/ 2 , (21) 

i.e., under typical coronal conditions for non-flaring loops (T < 10 MK, \ ~ 
—0.5), the DEM slope is near unity. If strong thermal conduction is included 
at the footpoints. then t he slope changes to +3/2 if not too close to T a 
l|van den Oord et alll997t) . The single-loop DEM sharply increases at T w T a 

(Fig- El- 
Loop expansion (a > 0) obviously steepens the DEM. Increased heat- 
ing at the loop footpoints (instead of uniform heating) makes the T range 
narrower and will also increase the slo pe of the DEM (see numerical cal- 
culations of various loop exa mples by ISchriiver k. Aschwandenl 120021 and 
lAschwanden fc Schriiverlho02l Fig. Eft). 

The DEM of flaring structures 

lAntiochosI l|l980|) discussed DEMs of solar flaring loops that cool by i) static 
conduction (without flows), ii) evaporative conduction (including flows), and 
iii) radiation. The inferred DEMs scale, in the above order, like 
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Qcond ^ T i Qevap &~ T , Qiad 7" ■ (^^) 

Since x « 0±0.5 in the range typically of interest for stellar flares (5—50 MK), 
all above DEMs are relatively flat (slope 1 ± 0.5). If multiple loops with equal 
slope but different peak T contribute, then the slope up to the first DEM 
peak can only become smaller. Non-constant loop cross sections have a very 
limited influence on the DEM slopes. 

Stellar flare observations are often not of sufficient quality to derive tem- 
perature and EM characteristic s for many differe nt time bins. An interesting 
diagnostic was presented by iMewe et alJ l|l997|) who calculated the time- 
integrated (average) DEM of a flare that decays quasi-statically. They find 

Q oc T 19 / 8 (23) 

up to a maximum T that is equal to the temperature at the start of the decay 
phase. 

Systems o f episodically flaring loops were computed semi-analytically by 
Cargill (1994), using analytic approximations for conductive and radiative 
decay phases of the flares. Here, the DEM is defined not by the internal 
loop structure but by the time evolution of a flaring plasma (assumed to be 
isothermal). Cargill argued that for radiative cooling, the (statistical) contri- 
bution of a flaring loop to the DEM is, to zeroth order, inversely proportional 
to the radiative decay time, which implies 

Q(T) oc T~ x+1 (24) 

up to a maximum T m , and a factor of T 1 / 2 less if subsonic draining of the 
cooling loop is allowed. Simulations with a uniform distribution of small flares 
within a limited energy range agree with these rough predictions, indicating a 
time-averaged DEM that is relatively flat below 10 6 K but steep (Q[T] oc T 4 ) 
up to a few MK, a range in which the cooling function drops rapidly. 

Let us next assume - in analogy to solar flares - that the occurrence 
rate of stellar flares is distributed in energy as a power law with an index 
a (dN/dE oc E~ a ) Then, an analytic expression can be derived for the 
time-averaged DEM of such a flare ensemble, i.e., a "flare-heated corona" 
l(Giidel et al.ll2003a| . We present a brief outline of the derivation. 

Observationall y, the flare peak tem perature is correla ted with the peak 
EM both in solar jFeldman et al.lll995|) and stellar flares (|Gudelll20o'^ : 

EM = aT b [cm -3 ] (25) 

with b w 4.3 ± 0.35 in the range of T = 10 — 100 MK (see Fig. |2"U1 further 
below) . The X-ray luminosity Lx of a plasma due to bremsstrahlung and line 
emission can be expressed as 



L x « EM A(T) = g EM 



(26) 



16 Manuel Giidel 



with <j) Ri 0.3 over the above temperature range (for broad-band X-ray losses). 

From (|25|1 and (|26[1 we obtain a relation between the flare peak tempera- 
ture T and its peak luminosity Lxo, 

j \ 1/(6-*) 



T a = . (27) 

V a 9 J 

We will investigate the general case in which r varies with the flare energy, 
namely 

r = roE? (28) 

(where /3 > is assumed, and To is a constant adjusted to the larger detected 
flares) in which case 



dN = dN dE = k > L -(»~^ ' ' ' 
dLxs) dE dL x .r Xfi 



k , < (29) 



where the constant k' — fer^ 1- "^ ^ /(I — (3) > as long as f3 < 1 (which 
can be reasonably assumed). Since we neglect the short rise time of the flare, 
our flare light curves are described by their exponential decay at t > 0, 

Lx(t) =L Xfi e- l ' T . (30) 

From hydrodynamic modeling, theory and observations, it is known that 
durin g the flare decays Tocn^, where n e is the plasma density ijReale et al.l 
119931) . The parameter ( is usually found between 0.5 and 2. 

Integrating the above equations for all emission contributions from the 
decaying and cooling flare plasma over the entire distribution of flares leads 
to two expressions valid for different regimes: 

f T 2 /^ for T < T m 

Q( T ) °C j T -(b_0)( Q _ 2 /3)/(l_/3)+2(,-0 for T>tZ ^ 31 ) 

Here, T m (a free parameter) is the temperature of the DEM peak. It is con- 
trolled by t he lower cutoff in t he power-law energy distribution of flares (for 
details, see lGiidel etal]l2003a|) - 



5.5 Discussion and summary on the temperature structure 

There is, as of yet, no clear explanation for the observed DEMs in active 
stars. The most notable feature is the often steep slope on the low-T side. The 
steepness is not easy to explain with standard static loop models, although 
strong footpoint heating may be a way out. Another explanation are magnetic 
loops with an expanding cross section from the base to the apex. The larger 
amount of plasma at high temperatures near the apex evidently steepens the 
DEM. Solar imaging does not prefer this type of loop, however. 
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Alternatively, the steep slope may be evidence for continual flaring; equa- 
tion l|31[l predicts slopes between 1 and 4, similar to what is often found in 
magnetically active stars. There is other evidence for this model to have some 
merit - we will encounter it again in the subsequent sections. 

It is also well established that more luminous stars (of a given size) reveal 
hotter coronae (Fig. UJ. Again, the cause of this relation is unclear. Possi- 
bly, increased magnetic activity leads to more numerous interactions between 
adjacent magnetic field structures. The heating efficiency thus increases. In 
particular, we expect a higher rate of large flares. The increased flare rate pro- 
duces higher X-ray luminosity because chromosphcric evaporation produces 
more EM; at the same time, th e plasma is heated to higher temperatures in 



5.6 Electron densities in stellar coronae 

Coronal electron densities control radiative losses from the coronal plasma; 
observationally, they can in principle also be used in conjunction with EMs 
to derive approximate coronal source volumes. The spectroscopic derivation 
of coronal densities is subtle, however. Two principal methods are available. 

Densities from Fe line ratios: The emissivities of many transitions of Fe 
ions in the EUV range are sensitive to densities in the range of interest to 
coronal research (jMewe et alJll9 85). The different density dependencies of 
different lines of the same Fe ion then also make their line-flux ratios, which 
(apart from blends) are easy to measure, useful diagnostics for the electron 
density. 

Densities from He-like line triplets: The He-like triplets of Cv, Nvi, 
Ovil, Neix, Mgxi, and Sixni provide another interesting density diagnostic 
for stellar coronae. Two examples are shown in Fig. (right). The spectra 
show, in order of increasing wavelength, the resonance, the intercombination, 
and the forbidden line of the O vil triplet. The corresponding transitions 
are depicted in the left panel of the Figure. The ratio between the fluxes 
in the forbidden line an d the intercombination line is sensitive to density 
( Gabriel & Jordan! Il96flft for the following reason: if the electron collision 
rate is sufficiently high, ions in the upper level of the forbidden transition, 
ls2s 3 Si, do not return to the ground level, Is 2 1 S'o, instead the ions are 
collisionally excited to the upper levels of the intercombination transitions, 
\s2p 3 -Pi, 2, from where they decay radiatively to the ground state. They thus 
enhance the flux in the intercombination line and weaken the flux in the 
forbidden line. The measured ratio R = f ji of the forbidden to the intercom- 
bination line flux can be written as 




R 



Rq 



f 



(32) 



1 + n e /N c 



where Rq is the limiting flux ratio at low densities and N c is the critical 
density at which R drops to Rq/2. The parameters Rq and N c are slightly 
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dependent on the electron temperature in the emitting source. A few use- 
ful parameters are collected in Tabled A systematic problem with He-like 
triplets is that the critical density N c increases with the formation tempera- 
ture of the ion, i.e., higher- Z ions measure only high densities at high T, while 
the lower-density analysis based on Cv, Nvi, Ovn, and Neix is applicable 
only to cool plasma. 

* co11 "'™ OVII triplets of Algol and Capella 




X / A 



Fig. 8. Left: Term diagram for transitions in He-like triplets. The resonance, inter- 
combination, and forbidden transitions are marked. The transition from 3 Si to 3 Pi 
re-distributes electrons from the upper level of the forbidden transition to the upper 
level of the intercombination transition, thus making the f/i line- flux ratio density 
sensitive. In the presence of a strong UV field, however, the same transition can be 
induced by radiation as well. Right: He- like triplet of O VII for Capella (black) and 
Algol (green). The resonance (r), intercombination (i), and forbidden (f) lines are 
marked. The f/i flux ratio of Algol is suppressed probably due to the strong UV 
radiation field of the primary B star (data from Chandra; both figures courtesy of 
J.-U. Ness). 



A review of the literature shows a rather unexpected segregation of coronal 
densities into two realms at different temperatures. The cool coronal plasma 
measured by C, N, and O lines in inactive stars is typically found at low, 
solar-like densities of order 10 9 cm -3 - 10 10 cm -3 . In active stars, the cooler 
components may show elevated densities up to several times 10 10 cm -3 , but 
it is the hotter plasma c omponent that apparently reveals extrem e values 
up to > 10 13 cm" 3 (e.g.. lDunree et alJll993 ISchriiver et alJll995lk A basic 
concern with these latter measurements is that most of the reported densities 
are only slightly above the low-density limits for the respective ratios, and 
upper limits have equally been reported, sometimes result ing in conflicting 
statements for different line ratios in the same spectrum l|Ness et al.l f2004; 
IPhillips et al]l200li lOsten et al.ll2003j) . Several authors have concluded that 
the extremely high densities found in some active stars are spurious and per- 
haps not representative of coronal features. The observational situation is 
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Table 2. Density-sensitive He-like triplets" 



Ion 


A(r,i,/) (A) 


Ro 


N c 


logrie range 6 


T range c (MK) 


Cv 


40.28/40.71/41.46 


11.4 


6 x 10 s 


7.7-10 


0.5-2 


Nvi 


28.79/29.07/29.53 


5.3 


5.3 x 10 9 


8.7-10.7 


0.7-3 


Ovu 


21.60/21.80/22.10 


3.74 


3.5 x 10 10 


9.5-11.5 


1.0-4.0 


Neix 


13.45/13.55/13.70 


3.08 


8.3 x 10 11 


11.0-13.0 


2.0-8.0 


Mg xi 


9.17/9.23/9.31 


2.66 d 


1.0 x 10 13 


12.0-14.0 


3.3-13 


Si xiii 


6.65/6.68/6.74 


2.33 d 


8.6 x 10 13 


13.0-15.0 


5.0-20 



"data derived from Porauct et al. (2001) at maximum formation temperature of ion 
6 range where R is within approximately [0.1,0.9] times Ro 
c range of 0.5-2 times maximum formation temperature of ion 
d for measurement with Chandra HETGS-MEG spectral resolution 



clearly unsatisfactory at the time of writing. The resolution of these contra- 
dictions requires a careful reconsideration of atomic physics issues. 

6 The structure of stellar coronae 

The magnetic structure of stellar coronae is one of the central topics in the 
stellar coronal research discipline. The extent and predominant locations of 
magnetic structures currently hold the key to our understanding of the inter- 
nal magnetic dynamo. All X-ray inferences of coronal magnetic structure in 
stars other than the Sun are so far indirect, while direct imaging, although 
at modest resolution, is available at radio wavelengths. 

6.1 Magnetic loop models 

Closed magnetic loops are the fundamental "building blocks" of the solar 
corona. When interpreting stellar coronae of any kind, we assume that this 
concept applies as well, although caution is in order. Even in the solar case, 
loops come in a wide variety of shapes and sizes (Fig. [SJ and appear to imply 
heating mec hanisms and heatin g locat ions that are poorly understood - see, 
for example. lAschwanden et all l|200fj|) . Nevertheless, simplified loop models 
offer an important starting point for coronal structure studies and possibly for 
coronal heating diagnostics. A short summary of some elementary properties 
fol lows. 

iRosner et alJ l)l978f) (RTV) have modeled hydrostatic loops with constant 
pressure (i.e., the loop height is smaller than the pressure scale height). They 
also assumed constant cross section, uniform heating, and absence of gravity, 
and found two scaling laws relating the loop semi-length L (in cm), the 
volumetric heating rate e (in erg cm _3 s~ 1 ), the pressure p (in dynes cm -2 ), 
and the loop apex temperature T a (in K), 
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Fig. 9. Left: Example of a solar coronal loop system observed by TRACE. Right: 
Flaring loop system (observation by TRACE at 171 A). Although these images 
show the emission from relatively cool coronal plasma, they illustrate the possible 
complexity of magnetic fields. 



T a = 1400(pi) 1/3 ; e = 9.8 x 10 4 / /6 L~ 5/6 . (33) 

ISerio et alJ lll98lf) extended these scaling laws to loops exceeding the pressure 
scale height s p , whereby, however, the limiting height at which the loops grow 
unstable is (2 — 3)s p : 

T a = 1400(p£) 1/3 e~°- 04i(2/sH+1/s *' ) ; e = 10y /6 £~ 5/6 e°- 5L(1/sff ~ 1/Sp) . 

^ ' / ^ (34) 

Here, sh is th e heat deposition scal e height. For loops with an area expansion 
factor r > 1. IVeseckv et alJ lll979l) found numerical solutions that approxi- 
mately follow the scaling laws jSchriiver et alJll989bh 

T a « U00r- o l (pL)^ 3 ; T a = 60r- 01 L 4 / 7 e 2 / 7 . (35) 

There are serious disagreements between some so/ar-loop observations and 
the RTV formalism so long as simplified quasi-static heating laws are as- 
sumed, the loops being more isothermal than predicted by the models. 
There is, however, only limited understanding of possible remedies, such 
as heating that is strongly concentrated at the loop footpoints, or dynam- 
ical processes in the loops (se e, for example, a summary of this debate in 
ISchriiver fc Aschwandenlbo02(l . 

6.2 Coronal structure from loop models 

When we interpret stellar coronal spectra, we assume, to first order, that some 
physical loop parameters map on our measured quantities, such as tempera- 
ture and EM (and possibly density), in a straightforward way. In the simplest 
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approach, we assume that the observed luminosity Lx is produced by an en- 
semble of identical coronal loops with characteristic half-length L, surface 
filling factor /, and an apex temperature T used for the entire loop; then, on 
using (|33|) and identifying Lx = eV, we obtain 



This relation can only hold if L is smaller than the pressure scale height. 
As an example, for an active solar analog (R — R@, Lx = 10 30 erg s -1 , 
T = 10 7 K) we obtain L ps 2x 10 11 / cm. The coronal volume is approximately 
V ps 8R 2 fL and the electron density n e = (EM/V) 1/2 . Further, L x « 
2x 10~ 23 EM erg s _1 . For the solar analog, thus, we find n e f pa 2.5x 10 9 cm -3 . 
The luminous, hot plasma component in magnetically active stars therefore 
seems to invariably require either very large, moderate-pressure loops with 
a large filling factor, or solar-sized high-pressure compact loops with a very 
small ( <, 1%) filling factor. 

While the above interpretational work identifies spectral-fit parameters 
such as T or EMs with parameters of theoretical loop models, a physically 
more appealing approach involves full hydrostatic models whose calculated 
emission spectra are directly fitted t o the observations. 

S uch studies dStern et alJIlQSfit iGiamnana et alJIlQQfit iMaggio fc Peresl 
Il997t IVentura et alJll998l) have found that the cooler component at « 1 — 
2 MK requires loops of small length (L but high pressure (p > p Q ), 

whereas the high-T component at ps 5 — 10 MK must be confined by very 
compact loops with extremely high base pressures (up to hundreds of dynes 
cm~ 2 ) and small (< 1%) filling factors. These parameters are suspiciously 
"flare-like" - the observed hot plasma is perhaps indeed related to multiple, 
very compact flaring regions. 

The most essential conclusion from these exercises is perhaps that, within 
the framework of such simplistic models, the loop heating rate required for 
magnetically active stars may exceed values for typical solar loops by orders 
of magnitude, pointing toward some enhanced heating process reminiscent 
of the energy deposition in flares. The compactness of the hot loops and the 
consequent high pressures also set these coronal structures clearly apart from 
any non-flaring solar coronal features. 

6.3 Coronal structure from densities 

Spectroscopically measured densities provide, in conjunction with EMs, im- 
portant estimates of emitting volumes. If the trend suggested from density- 
sensitive line flux ratios holds, namely that for increasing temperature, the 
pressures become progressively higher, then progressively smaller volumes 
are a consequence. The volume required for a luminosity of 10 30 erg s _1 at 
10 MK is V = L x /(2 x 10" 23 n 2 ) cm" 3 = 5 x 10 26 cm" 3 (where the coefficient 




[cm]. 



(36) 
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2 x 10 -23 is from (J5J, appropriate for T = 10 MK), corresponding to a layer 
of only 80 m height around a solar-like star, or 8 km for a filling factor of 
only 1%! Such scales are much smaller than chromospheric scale height and 
therefore problematic. Still smaller filling factors must be assumed for a star 
of this kind. Similarly, from the RTV loop scaling law l|33|) - if applicable 
-, a loop height h = 2L/ir — 8.5 x 10 5 T 2 /n e ~ 80 km is found, again an 
unreasonably small size. 

The confinement of plasma at such high densities in compact sources 
would also require coronal magnetic field strengths of order B > (lQiniekT) 1 / 2 
w 1 kG, i.e., field strengths like those very close and just above (sun-)spots. 
In that case, the typical magnetic dissipation time is only a few seconds for 
n e w 10 13 cm" 3 if the energy is derived from the same magnetic fields, sug- 
gesting that the small, bright loops light up only briefly. In other words, the 
stellar corona would be made up of numerous ephemer al loop sources that 
canno t be treated as being in a quasi-static equilibrium l|van den Oord et all 
119971) . 

6.4 X-ray coronal imaging: Overview 

X-ray images of stellar coronae have been derived from eclipses in binaries, or 
from rotational modulation in rapidly rotating stars. We keep in mind that 
any indirect imaging of this kind is highly biased by observational constraints 
(e.g., the volume that is subject to eclipses or self-eclipses, or the accessible 
temperature range) and by the amount and density of plasma trapped in the 
magnetic fields. X-ray imaging captures strongly emitting plasma, not the 
entire magnetic-field structure. 

The "image" to be reconstructed consists of volume elements at coordi- 
nates (x,y,z) with optically thin fluxes f(x,y,z) assumed to be constant in 
time. In the special case of negligible stellar rotation during the observation, 
the problem can be reduced to a 2-D projection onto the plane of the sky, at 
the cost of positional information along the line of sight (Fig. 1101) . In general, 
one thus seeks the geometric brightness distribution f(x,y,z) = /y/. (i,j,k 
being the discrete number indices of the volume elements) from a binned, ob- 
served light curve F s = F(t s ) that undergoes a modulation due to an eclipse 
or due to rotation. 

6.5 Active-region modeling 

In the most basic approach, the emitting X-ray or radio corona can be mod- 
eled by making use of a small number of simple, elementary building blocks 
that are essentially described by their size, their brightness, and their loca- 
tion. This approach is the 3-D equivalent to standard surface spot modeling. 
Preferred building block shapes are radially directed, uniformly bright, opti- 
cally thin, radially truncated spherical cones with their apexes at the stellar 
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Fig. 10. S ketch showing the g eometry of an eclipsing binary (in this case, a CrB, 
figure from lGudel et al.l2003lJ) . The large circles illustrate the limbs of the eclipsing 
star that moves from left to right in front of the eclipsed star (shown in yellow). 
The limbs projected at different times during ingress and egress define a distorted 
2-D array (x, y) of pixels (an example of a pixel is shown in gray). 



center. Free parameters are their opening angles, their heights above the stel- 
lar surface, their radiances, and their central latitudes and longitudes. These 
parameters are then varied until the model fits the observed light curve. 

A minimum solutio n can be found for a rotationally modulated star 
dOudel fc Schmittlll995lb If a rotationally modulated feature is invisible dur- 
ing a phase interval tp of the stellar rotation, then all sources contributing to 
this feature must be confined to within a maximum volume, V maxi given by 

Vmax = ± _ (27r-y)(l+sin 2 z) + 2cot(x/2) 
Rl 3 6sini 3tani 

where tanx = tan(</?/2)cosi, sm(i/j/2) = sin(c/?/2)sim with < ip/2 < tt/2, 
and x an d f/2 lie in the same quadrant (i is the stellar inclination, < i < 
7r/2). Together with the modulated fraction of the luminosity, lower limits to 
average electron densities in the modulated region follow directly. 
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6.6 Maximum-entropy image reconstruction 

Maximum entropy methods (MEM) are applicable both to rotationally mod- 
ulated light curves and to eclipse observations. The standard MEM selects 
among all images fijk (defined in units of counts per volume element) that 
are compatible with the observation, the one that minimizes the Kullback 
contrast ("relative entropy") 

* = (38) 

i,j,k l% i k 

with respect to an a priori image f^ k , which is usually unity inside the allowed 
area or volume and vanishes where no brightness is admitted. Minimizing K 
thus introduces the least possible information while being compatible with 
the observation. The contrast K is minimum if fy^ is proportional to /°- fc 
and thus flat inside the field of view, and it is maximum if the whole flux is 
concentrated in a single pixel (i,j,k). The compatibility with the observed 
count light curve is measured by x 2 , 

X 2 =E^=^ (39) 

where F s and F* are, respectively, the observed number of counts and the 
number of counts predicted from f^k and the eclipse geometry. Poisson statis- 
tics usually requires more than 15 counts per bin. Finally, normalization is 
enforced by means of the constraint 



* = ~" (40) 

where f tot is the sum of all fluxes in the model. 
The final algorithm minimizes the cost function 

C = x 2 +tK + rjN. (41) 

The trade-off between the compatibility with the observation, normalization, 
and unbiasedness is determined by the Lagrange multipliers £ and r\ such 
that the reduced y 2 is and normalization holds within a few percent. 



6.7 Lucy/Withbroe image reconstruction 

This method ( after lLucvlll974l and lWithbroelli975l Uiterativelv adjusts fluxes 
in a given set of volume elements based on the mismatch between the model 
and the observed light curves in all time bins to which the volume element 



Stellar Coronae 25 



contributes. At any given time t s during the eclipse, the observed flux F s is 
the sum of the fluxes fijk from all volume elements that are unocculted: 



where m s is the "occultation matrix" for the time t s : it puts, for any given 
time t s , a weight of unity to all visible volume elements and zero to all 
invisible elements (and intermediate values for partially occulted elements). 
Since F s is given, one needs to solve (|42|) for the flux distribution, which is 
done iteratively as follows: 



where F (t s ) and F^(t s ) are, respectively, the observed flux and the model 
flux (or counts) in the bin at time t s , both for the iteration step n. Initially, 
a plausible, smooth distribution of flux is assumed, e.g., constant brightness, 
or some r~ p radial dependence. 

6.8 Backprojection and Clean image reconstruction 

If rotation can be neglected during an eclipse, for example in long-period de- 
tached binaries, then the limb of the eclipsing star is projected at regular time 
intervals onto the plane of the sky and therefore onto a specific part of th e 
eclipsed corona, first during ingress, later during egress llOudel et al]l2003bh . 
The two limb sets define a 2-D grid of distorted, curved pixels (Tig. llOfl . The 
brightness decrement during ingress or, respectively, the brightness increment 
during egress within a time step [t s , i s +i] originates from within a region con- 
fined by the two respective limb projections at t s and t s+ \. Ingress and egress 
thus each define a 1-D image by backprojection from the light curve gradients 
onto the plane of the sky. The relevant reconstruction problem from multiple 
geometric projections is known in tomography. The limiting case of only two 
independent projections can be augmented by a CLEAN step, as follows. 
The pixel with the largest sum of projected fluxes from ingress and egress is 
assumed to represent the location of a real source. A fraction, g < 1, of this 
source flux is then subtracted from the two projections and saved on a clean 
map, and the process is iterated until all flux is transferred onto the latter. 

6.9 X-ray coronal structure inferred from eclipses 
Extent of eclipsed coronal features 

Some shallow X-ray eclipses in tidally interacting binary systems of the RS 
CVn, Algol, or BY Dra type have provided important information on ex- 
tended coronal structure. For example, Walter et al. (1983) concluded that 




(42) 




(43) 



S 
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the coronae in the AR Lac binary components are bi-modal in size, consist- 
ing of compact, high-pressure (i.e., 50-100 dynes cm~ 2 ) active regions with 
a scale height < i?„ , while the subgiant K star is additionally surrounded by 
an extended (2.7i?») low-pressure corona. Further, a hot component pervad- 
ing the entire binary system was im plied from the abs ence of an eclipse in 
the hard ME detector on EXOSAT ijwhite et al.lll990h . and simi lar conclu- 
sions have been drawn from detailed light-curve inversion analysis l|Pres et alJ 

X-ray dips or any periodic modulation have often been absent in X-ray 
observations of binaries. This again has b een taken as evid ence of a very 
extended (about li?* in the case of Algol, IWhite et al.lll98fil) X -ray corona 
unless more compact structures sit at high latitudes where they remain un- 
eclipsed. 

Among wide, non-interacting eclipsing stars, a CrB provides a particu- 
larly well-suited example because its X-ray active, young solar analog (G5 V) 
is totally eclipsed every 17 days by the optical primary, an AO V star that 
is perfectly X-ray dark. Other parameters are ideal as well, such as the non- 
central eclipse, the eclipse time-scale of a few hours, and the relatively slow 
rotati on period of the seco ndary. Image reconstructions from eclipse observa- 
tions l|Gudel et alJl2003b|) reveal patches of active regions across the face of 
the G star; not much material is found significantly beyond its limb (Fig. 1120. 
The structures tend to be of modest size (»5x 10 9 cm), with large, X-ray 
faint areas in between, although the star's luminosity exceeds that of the 
active Sun by a factor of «30. These observations imply moderately high 
densities in the emitting active regions, reaching a few 10 10 cm -3 in the 
brightest active regions. 

6.10 Eclipsed X-ray flares 

Eclipses of flares have contributed very valuable information on densities and 
the geometric size of flaring structures. Only few reports are available, among 
th em the following : 

IChoi fe Dot anil (1998) described a full eclipse of an X-ray flare in progress 
in the contact binary system VW Cep. During a narrow dip in the flare 
decay, the X-ray flux returned essentially to the pre-flare level. Geometric 
considerations then placed the flare near one of the poles of the primary 
star, with a size scale of order 5.5 x 10 10 cm or somewhat smaller than 
the secondary star. The authors consequently inferred an electron density 
of 5 x 10 10 cm~ 3 . A polar loc ation was also advocated for a flare on Algol 
observed across an eclipse bv ISchmitt fc Favatal l|l999h . The flare emission 
was again eclipsed completely, and judged from the known system geometry, 
the flare was located above one of the poles, with a maximum source height 
of no more than approximately 0.5i?*, implying a minimum electron density 
of 9.4 x 10 10 cm -3 if the volume filling factor was unity. A more moderate 
flare was observed during an eclipse in the Algol system by ISchmitt et alJ 
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Fig. 11. Two examples of eclipses and the corresponding coronal imag e reconstruc- 
tions. From top to bottom: Light curve of the YY Gem system (from iGiidel et alJ 
l2001al. observation with XM M-Newton EPIC); light curve of the AR Lac system 
( after ISiarkowski et al.lll9 96. observation with ASCA SIS); reconstructed image of 
the coronal structure of, respectively, YY Gem (at phase 0.375) and AR Lac (at 
quadrature). The latter figure shows a solution with intrabinary emission. (The 
light curve of AR Lac is phase-folded; the actual observation started around phase 
0; data and image for AR Lac courtesy of M. Siarkowski.) 



(2003). In this case, the image reconstruction required an equatorial location, 
with a compact flare source of height h « O.li?*. Most of the source volume 
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Fig. 12. Light curve and image reconstruction of the A+G binary a CrB. The 
left panel shows the light curve from observations with XMM-Newton, the right 
panel illustrates the reconstructed X-ray brightness distribution on the G star (after 
iGiidel et al.ll2003rl . 



exceeded densities of 10 11 cm" 3 , with the highest values at ~ 2 x 10 11 cm -3 . 
Because the quiescent flux level was attained throughout the flare eclipse, 
the authors argued that its source, in turn, must be concentrated near the 
polar region with a modest filling factor of / < 0.1 and electron densities of 
w 3 x 10 10 cm~ 3 . 



6.11 Radio Very Long Baseline Interferometry 

By interferometrically combining radio telescopes over large distances, an- 
gular resolutions of as little as one milliarcsecond can be achieved in the 
microwave range. 

VLBI techniques have been very demanding for single late-type dwarf 
stars, owing both to low flux levels and small coronal sizes. Some observations 
with mas angular resolutions show unresolved quiescent or flari ng sources, 
thus constraining the brightness temperature to If, > 10 10 K (e.g.. lBenz et alJ 
whereas others show evidence for extended coronac with coronal sizes 
up to several times the stellar size. 

The dMe star UV Cet was found to be surrounded by a pair of giant syn- 
chrotron lobes, with sizes up to 2.4 x 10 10 cm and a separation of 4—5 stellar 
radii along the putative rotation axis of the star, suggesting very extended 
magnetic structures above the magnetic poles fFig. H3b ,'l. perhaps arranged in 
a global dipol e as sketched in Fi g. ITBb . I discuss this observation is some de- 
tail, following nz et al.l l)l998(l . to demonstrate the procedures with which 
we can characterize the magnetic field structure. Throughout, we assume, 
as detailed further above, that the emission is optically thin gyrosynchrotron 
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emission from a power-law population of accelerated electrons with a number 
density distribution in energy e 

n(e) = N(5 - [cm^erg" 1 ] (44) 

where e = (7 — l)m e c 2 is the kinetic particle energy, 7 is the Lorentz factor, 
and S > 1 has been assumed so that N is the total non-thermal electron 
number density above eo- 

We first need an expression for the emissivity: For isotropic pitch an- 
gle electron distributions according to l|44ll with 2 <S 5 7, for harmonics 
10 < s < 100, and for the x-mode <|Dulklll985[) 

/ 1.22-0.905 

T) U « 3.3 x 10- 24 i?7V lO- o - 525 (sin0)-°- 43 + - 655 (^) (45) 

(in erg s _1 cm -3 Hz" 1 sterad" 1 ). Here, 9 denotes the emission angle to 
the magnetic field, and v c = eB/m e c is the electron gyrofrequency. The 
coefficient of this expression is applicable if a power-law cutoff has been set 
at 10 keV. This cut-off itself is of little relevance for gyrosynchrotron emission 
since electrons radiate little at such low energies. 

We assume = ir/2 and 6 « 2.5 (the latter from modeling of M dwarf 
spectra, Sect. l4~3*)l . Adding o- and x-mode (which have similar emissivities), 
the intensity is 

/ w 1.8t) v D (46) 

which is to be compared with the observed intensity of 3.9 x 10 -8 erg -1 
cm~ 2 Hz -1 sterad -1 . Here, the fitted diameter of one of the radio blobs 
(0.3 milliarcsec or D = 1.2 x 10 10 cm) was used. Combining l|45(l - H46|l . we 
find 

B«2x 10 5 iV~ 1/2 . (47) 

This relation is shown in Fie . 113b together with the condition that the particle 
pressure in the magnetic loops is less than the magnetic pressure, 

ATe-<g. (48) 

We find a lower limit to B of approximately 15 G and an upper limit to AT 
of about 2 x 10 8 cm -3 . 

At least part of the observed emission was slowly decaying, with a decay 
timescale of r = 6650 s (perhaps from a flare that filled the magnetospheric 
volume with electrons). Given the large size of the structure, the dominant 
energy decay process is likely to be due to synchrotron radiation loss (|T2)l . 
The average frequency of synchrotron emission can be expressed as 

v = 1.3 x 10 6 B7 2 [Hz] (49) 

l)Melroselll98flfl . Using the observing frequency (8.4 GHz in this case) for P, 
we find the decay time 
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log electron density [cm 3 ] 

Fig. 13. (a) VLBA image of the dMe star UV Cet; the two radio lobes are separated 
by about 1.4 mas, while the best angular resolution reaches 0.7 mas. The straight 
line shows the orientation of the putative rotation axis, assumed to be parallel to 
the axis of the orbit of UV Cet around the nearby Gl 65 A. The small circle gives 
the photospheric diameter to size, although the precise position is unknown (after 
Benz et al 1998). (b) Estimates of the magnetic field and the nonthermal electron 
density for the same observation (see text for details). 



r = 2 ~ 8 x 10 6 B- 3/2 [s] (50) 

7 

and therefore B = 113 G. This defines an upper limit to the magnetic field 
strength because we have ignored other energy losses that might be present. 
The upper limit is also drawn in Fig. 113b. We have thus confined the magnetic 
field strength in the source to 15-130 G. 

VLBA imaging and polarimetry of Algol reveals a similar picture with 
two oppositely polarized radio lobes separated along a lin e perpendicular to 
the orbital plane by more than the diameter of the K star ijMutel et al.l [l998. 
Fig. O an d Fig. 115b). Large-scale polarization structure further supports 
model s that assume globally organized magnetic fields around active binary 
stars l|Beaslev fc GudefeOOOj) . 

An important, early VLBI result for RS CVn and Algol-like binaries is 
evidence for a compact core plus extended halo r adio structure of a total 
size that is comparable to the binary system size (jMutel et alJll985)l . The 
basic idea here is the following: During quiescence, an optically thin, very 
large magnetosphere is filled with a power-law distribution of electrons. The 
emission is essentially optically thin, with a flat radio spectrum. During an 
outburst, an active region injects a larger electron population into magnetic 
loops. The compact source (the "core"), unresolved by VLBI, becomes op- 
tically thick due to synchrotron self-absorption, and the radio spectral in- 
dex consequently becomes positive, and the brightness temperature is thus 
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Fig. 14. VLBA observation of Algol at 8.4 GHz, resolving two lobes around the 
binary. The most likely configuration of the binary components is also drawn (after 
iMutel et al.lll99Sl) . 



equal to the effective electron temperature. IMutel et al.l l|1985|) measured a 
few times 10 10 K for outbursts on RS CVn-type binaries. Once the electron 
injection has ceased, the lifetime of the energetic electron population is de- 
termined essentially by the synchrotron loss time. The final phase consists 
of magnetic-loop expansion owing to buoyancy, expanding the source until it 
merges with the pre- flare "halo" component. At the same time, the source 
becomes progressively more optically thin, developing a negative radio spec- 
tral index and mild circular polarization. Given the large size, the strength of 
the extended magnetic field must be rather moderate, of order 10 G. To emit 
radiation at frequency v in the microwave range, high Lorentz factors are re- 
quired. Because the spectral power for synchrotron emission is predominantly 
emitted at a harmonic of order 7 3 of the relativistic gyrofrequency 



we require 7 6 x 10 _4 (^/_B) 1 / 2 . For an observing frequency of 5 GHz thus, 
7 « 10. A model calculation of the spectral development of the microwave 
radiation was outlined in Sect. 14.31 

6.12 Radio magnetospheric models 

VLBI observations of RS CVn and Algol binaries, T Tau stars, and magnetic 
Bp/Ap stars have shown some perplexing structures with sizes at least as 
large as the binary system, with polarization properties that suggest that the 
magnetic fields are globally ordered. This has led to a series of large-scale 
magnetic models for such systems. They have in common a global, dipole- 
like structure somewhat resembling the Earth's Van Allen belts (Fig. 115k). 



eB 



= 2.8 x 10 6 — , 

7 



(51) 



m e jc 
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Fig. 15. (a) Equatorial model for the magnetosphere of the young B star SI in p 
Oph iAndre et al ] ll98Sh . (b) Sketch for radio emissi on from a global d ipole consis- 
tent with the VLBI observation shown in Fig. 113b llMutel et alJll998l : reproduced 
with permission of the A AS.) 



Stellar winds escaping along magnetic fields draw the field lines into a cur- 
rent sheet configuration in the equatorial plane. Particles can be acceler- 
ated in that region. They subsequently travel back to become trapped in the 
dipolar-like, equatorial magnetospheric cavit y. Variants of this radiatio n-belt 
model, partly based on theoretical wor k of lHavnes fc Goert d l|l984lh have 
been applied to RS CVn bina ries (e.g.. | M oms et al.lll99dt) 7in an o ptically 
thick version to Bp/Ap stars ijDrake et alJll987t iLinskv et alJll992|) and in 
an optically thin version to a young B star ijAndre et alJll988|) . 

Polarization observations of RS CVn binaries support these models. The 
polarization degree seems to be anticorrelated with the luminosity for any 
given system, but the sense of polarization changes between lower and higher 
frequencies. For the entire binary sample, the polarization degree is inversely 
correlated with the stellar inclination angle such that low- inclination ( "pole- 
on" ) systems sh ow the strongest polarization degrees, as ex pected from such 
global systems l|Mutel et al.lll987Ul998tlMorris et al.lll990|) . 

Flares still appear to originate in compact sources within such magne- 
tospheres, probably close to the star. The "core plus halo" model then cor- 
rectly describes the radio spectral properties - the halo corresponds to the 
extended magnetosphere. This is also supported by magnetic fields inferred to 
be stronger (80—200 G) in the flaring core and weaker (10—30 G ) in the halo 
( ITmana et al Jl 993lll 99<*lMiit,fil et a .lll99fiHTrigilio et a.ll200ll> . Because the 
optical depth is frequency-dependent, the source is small at high frequencies 
(with a size above 10 GHz) and large at smal l frequencies (with a size 

comparable to the bin ary system size at 1.4 GHz; iKlein fc Chiuderi-Dragol 
ll987tljones et al.ll 995). This effect explains the relatively flat, optically thick 
radio spectra seen during flares. 
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6.13 Extended or compact coronae? 

As the previous discussions imply, we are confronted with mixed evidence for 
predominantly extended (source height > i?*) and predominantly compact 
(<C i?*) coronal structures or a mixture thereof, results that variably come 
from radio or from X-ray astronomical observations. There does not seem to 
be unequivocal agreement on the type of structure that generally prevails. 
Several trends can be recognized, however, as summarized below. 

Compact coronal structure. Steep (portions of the) ingress and egress light 
curves or prominent rotational modulation unambiguously argue in favor of 
short scale lengths perpendicular to the line of sight. Common to all are 
relatively high inferred densities (« 10 10 cm -3 ). The pressures of such active 
regions may exceed pressures of non-flaring solar active regions by up to 
two orders of magnitude. Spectroscopic observations of high densities and 
loop modeling add further evidence for the presence of some rather compact 
sources. Flare modeling also provides modest sizes, often of order 0.1-li?*, 
for the involved magnetic loops (Sect.lHJ. 

Extended structure. Here, the arguments are less direct and are usually 
based on the absence of deep eclipses, or very shallow ingress and egress 
curves. Caution is in order in cases where the sources may be located near 
one of the polar regions; in those cases, eclipses and rotational modulation 
may also be absent regardless of the source size. Complementary information 
is available from flare analysis (see Sect.[8J) that in some cases does suggest 
quite large loops. The caveat here is that simple single-loop models may not 
apply to such flares. Clear evidence is available from radio interferometry 
that proves the presence of large-scale, globally ordered magnetic fields. The 
existence of prominent extended, closed magnetic fields on scales > i?» is 
therefore also beyond doubt for several active stars. 

The most likely answer to the question on coronal structure size is there- 
fore an equivocal one: Coronal magnetic structures follow a size distribution 
from very compact to extended ( <; i?») with various characteristic densities, 
temperatures, non-thermal electron densities, and surface locations. Quite 
different loop systems may be responsible when measuring cool X-ray lines 
or hot lines, or when observing in the microwave range. This is no differ- 
ent from what we see on the Sun even though various features observed on 
magnetically active stars stretch the comparison perhaps rather too far for 
comfort. 



7 Stellar radio flares 

Flares arise as a consequence of a sudden energy release and relaxation pro- 
cess of the magnetic field in solar and stellar coronae. Present-day models 
assume that the energy is accumulated and stored in non-potential magnetic 
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fields prior to an instability that most likely implies reconnection of neigh- 
boring antiparallel magnetic fields. The energy is brought into the corona 
by turbulent footpoint motions that tangle the field lines at larger heights. 
The explosive energy release becomes measurable across the electromagnetic 
spectrum and, in the solar high-energy particles in interplanetary 

space as well. 

Flares are ubiquitous among coronal stars of all types, with very few 
exceptions. They have, of course, prominently figured in solar studies, and it is 
once again solar physics that has paved the way to the interpretation of stellar 
flares, even if not all features are fully understood yet. The complexity that 
flares reveal to the solar astronomer is inaccessible in stellar flares, especially 
in the absence of spatially resolved observations. Simplified concepts, perhaps 
tested for solar examples, must suffice. The following sections summarize the 
"stellar astronomer's way" of looking at flares. 

7.1 Incoherent radio bursts 

Active stars reveal two principal flare types at radio wavelengths, similar to 
the solar case: 

Incoherent flares evolve on time scales of minutes to hours, they show 
broad-band spectra and moderate degrees of polarization. These bursts are 
fully equivalent to solar microwave bursts. Like the latter, they show evidence 
for the presence of mildly relativistic electrons. The emission is therefore 
interpreted as gyrosynchrotron radiation from coronal magnetic fields. Many 
flares on single F/G/K s tars are of thi s type, as are almost all r adio flares on 
M dw arfs above 5 GHz ((BastiarJll99f^ . or on RS CVn binaries l(Mutel et alJ 

7.2 Coherent radio bursts 

The stellar equivalents to coherent solar radio bursts (showing high brightness 
temperature, short durations, perhaps small bandwidth, and perhaps high 
polarization degree) have been observed already in the early days of radio 
astronomy, given their sometimes extremely high fluxes. Like for the Sun, 
they come in a bewilderin g variety which has made a clear identification of 
the ultimate cause difficult flBastian et alJl998j) . Coherent bursts are frequent 
on late-t ype main-sequence stars, b ut have also been reported from RS CVn 
binaries ((White fc Franciosinilll995h . 

Such bursts carry important information in high-time resolution light 
curves. Radio "spike" rise times as short as 5—20 ms have been reported, 
requiring, from the light-time argument, source sizes of R < cAt ~ 1500- 
6000 km. With 

_ 2kT b v 2 ttR 2 , . 
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and measured fluxes up to 1 mJy for sources at a distance of a few parsec, 
we derive brightness temperatures of order T h sa 10 16 K, a clear proof of the 
presence of a coherent mechanism llLang et al.lll983l lLang fc Willsonlll98fit 
iGiidel et al.lll989at iBastian et al.lll990|) . 

7.3 Radio dynamic spectra 

The standard means to study solar coherent bursts are "dynamic spectra" 
(flux vs. frequency and time). If the elementary frequencies relevant for the 
coherent emission process, v p and v Cl evolve in the source, or if the radiating 
source itself travels across density or magnetic field gr adients, the emissio n 
leaves characteristic traces on the dynamic spectrum ijBastian et alJ ll998). 
The study of drifts, decay times, harmonic structures etc. may then help 
identify the emission process and thus infer magnetic field strengths, electron 
densities, electron energies and beam velocities etc. 

Applying the same technology to stellar observations has turned out to 
be extremely challenging, and only few successful dynamic spectra have been 
recorded to date. A rich phenomenology has been uncovered, including: a) 
short, highly polarized bursts with structures as narrow as Avjv = 0.2% 
suggesting either plasma emission from a sour ce of size ~ 3 x 10 8 cm, or a 

250 G JBastian fc Bookbinderlll987t 



cyclotron maser in 


magnetic fields 


) of ~ 


iGiidel et alJ 


1989a; 


; Bastian et alJ 


1990) 



with positive drift rates of 2 MHz s _1 around 20 cm wavelength, taken as 
evidence for a disturbance propaga ting "downward" in the plasma emission 
interpretation l|.Iackson et alll987lk and c) in solar terminology, rapid broad- 
band pulsations, "sudden (flu x) reductions" , and positive and negative drift 
rates of 250-1000 MHz s" 1 (|Bastian et al.lll990tlAb"ada-Simon et alJll99-d 

\\m% 

Spectral bandwidths as 1% of the emitting frequency were found for some 
bursts. If we conservatively assume a magnetic scale height of Lb = 1-R* (as- 
suming that the emission is gyromagnetic), the source size can be estimated 
to be 

r w —L B (53) 
v 

which is a few 1000 km for Lb = -R* of an M dwarf, aga in implying very 
high Th, compatible with the above light ti me argument ijLang fc Willsonl 
ll988UGudel et alJll989aHBastian et alJll990l) . 

Recent developments have permitted recording quite broad regions of the 
stellar burst spectrum, allowi n g for a much better characterization of the 
burst types. lOsten fc Bast ian (2006]) obtained burst spectra from AD Leo 
recorded in the 1120-1620 MHz range with a time resolution of 10 ms and a 
spectral resolution of 0.78 MHz. They find significant frequency-time drifts 
in short subbursts (Fig. 1170 . The main characteristics of these structures 
are (i) durations of about 30 ms, (ii) high polarization (> 90%), frequency 
bandwidths of Avjv « 5%, and inverse drift rates (time interval per drift 
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Fig. 16. Radio dynamic spectra of M dwarf flares. Upper three panels show a flare 
on AD Leo, recorded with the Arecibo (top), Effelsberg (mid dle) and Jodrell Ban k 
(bottom) telescopes in different wavelength ranges (see also iGiidel et al.l |f989a). 
Bottom three panels show fla res on AD Leo (top and middle) and YZ CMi (bottom) 
observed at Arecibo ( after IfBastian et alJll990L reproduced with permission of the 
AAS). 



in frequency) that are symmetric around zero, with a characteristic width 
corresponding to a frequency drift of 2.2 GHz s _1 . These characteristics are 
very close to those of narrow-band solar decimetric "spike bursts" but set 



Stellar Coronae 37 



them apart from type III bursts associated with electron beams in the solar 
corona (type III bursts have considerably longer durations, occupy broader 
frequency intervals, and are at best moderately polarized). 




131.0131.2131.4131.6131.8132.0132.2132.4 
Time from start (s) 

22:35:00 2E:35:20 23:35:40 23:3S:00 22:34:20 22:36:iO 
Slnrl Time (13-JUII-03 23:11:57) 

Fig. 17. Dynamic spectra of radio bursts observed in AD Leo dOsten fc Bastianl 
2006). The left figure shows the entire event, the right figure shows an extract 
revealing frequency drifts and spectral structure. 



There is one significant difference between the stellar and solar bursts: the 
duration of the solar examples is ab out three times shorter. We can relate 
this t o Coulomb collisional damping ijGiidel fe BendH99(j lOsten fe Bastiar] 
2006) which is determined by the electron-ion collision time r e i, from which 
we find 

T = 8100^ z r 2 / 3 . (54) 

We now identify r with the observed characteristic event duration, and 
with the radio fr equency in MHz. For solar events at 1.5 GHz, T w 5 MK 
( Gudel fe Ben3ll99d). but for the longer events on AD Leo, T « 13 MK 
( Osten fc Bastianll2006lh which can be e xplained by the higher coronal tem- 
perature of AD Leo l|Giidel et al.ll2003al) . Conversely, the higher temperature 
of AD Leo's corona increases the gyroresonance absorption coefficient for 
underlying electron cyclotron maser emission significantly: 

K gr cx T s_1 (55) 

where s = 2, 3 is the harmonic of the absorbing layer, making the escape of 
the preferred harmonics of the maser difficult. This is not the case for plasma 
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radiation for which free- free absorption is relevant: 

Kg oc T~ 3/2 , (56) 

implying that the higher temperature of AD Leo's corona clearly favors a 
plasma emission process as long as tu p > fl c . The latter condition is easily met 
as it requires B < 500 G. This is a nice example where a stellar observation 
helps identify the relevant emission mechanism that remains equivocal under 
solar conditions. 



8 Stellar X-ray flares 
8.1 Cooling physics 

Flares cool through radiative, conductive, and possibly also volume expansion 
processes. We define the flare decay phase as the episode when the net energy 
loss by cooling exceeds the energy gain by heating, and the total thermal 
energy of the flare plasma decreases. The thermal energy decay time scale 
T t h is defined as 

n h = | (57) 

where E rs 3n e kT (for n e w nn) is the total thermal energy density in the 
flaring plasma of electron density n e and temperature T, and E is the volu- 
metric cooling loss rate (in erg cm -3 s _1 ). For conduction across temperature 
gradients in parallel magnetic fields, the mean loss rate per unit volume is 

^ = ~ L ^-^—^ (58) 

where s is the coordinate along the field lines, and the term kqT 5 ^ 2 dT/ds 
is the conductive flux in the approximation of ISpitzerl l)l962j) . to be eval- 
uated near the loop footpoint where T drops below 10 6 K, with kq w 

9 x 10~ 7 erg cm _1 s _1 K -7 / 2 . Equations l|57l I58fl define the conductive time 
scale r t h = t c . The second equation in l|58|l should be used only as an ap- 
proximation for non-radiating loops with a constant cross section down to 
the loss region and with uniform heating (or for time-depende nt cooling of a 
constant-pressure loop with out heating; for the factor of 4/7, see lDowdv et al.l 
I1985HKqpp fc Polet,toll99l . We have used L for the characteristic dimension 
of the source along the magnetic field lines, for example the half-length of 
a magnetic loop. Strictly speaking, energy is not lost by conduction but is 
redistributed within the source; however, we consider energy lost when it is 
conducted to a region that is below X-ray emitting temperatures, e.g., the 
transition region/chromosphere at the magnetic loop footpoints. 

Radiative losses are by bremsstrahlung (dominant for T J> 20 MK), 2- 
photon continuum, bound-free, and line radiation. We note that the plasma 
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composition in terms of element abundances can modify the cooling function 
A(T), but the correction is of minor importance because stellar flares are 
usually rather hot. At relevant temperatures, the dominant radiative losses 
are by bremsstrahlung, which is little sensitive to modifications of the heavy- 
element abundances. The energy loss rate is 

E r = n e n„A(T) (59) 

For T > 20 MK, A(T) = ApT 1 sa iq-24.66 T i/4 erg cm 3 g -i ( after 
Ivan den Oord & Mewell989l and lMewe et al.ll985^ . Equations %?7\Wb define 
the radiative time scale r t h = T r . 



8.2 Interpretation of the decay time 

Equations J57J, JSHJ, and lf5^|) describe the decay of the thermal energy, 
which in flare plasma is primarily due to the decay of temperature (with a 
time scale tt) and density. In contrast, the observed light curve decays (with 
a time scale Td for the luminosity) primarily due to the decreasing EM and, 
to a lesser extent, due to the decrease of A(T) with decreasing temperature 
above w 15 MK. From the energy equation, the thermal energy decay time 
scale Tth is found to be 

-=fl-?) J - + ^ L (60) 
r th V 2J t t 2r d v ' 

where the right-hand si de is usually known from the observations (see 
Ivan den Oord et alJll988l for a deriv ation). The decay time scale of the EM 
then follows as 1/tem = 1/rw — ^/rr. lPan et al.l l|l997h derived somewhat dif- 
ferent coefficients in (IfiOJ) for the assumption of constant volume or constant 
mass, including the enthalpy flux. In the absence of measurements of tt, it 
is often assumed that r t h = Td although this is an inaccurate approximation. 

In (|60(l . T t h is usually set to be r r or t c or, if both loss terms are signif- 
icant, (r" 1 + t" 1 ) -1 , taken at the beginning of the flare decay (note again 
that a simple identification of r r with Td is not accurate). If radiative losses 
dominate, the density immediately follows from l|57l I59|l 

T th « -TT-r (61) 

n e A(T) 

and the characteristic size scale £ of the flaring plasma or the flare-loop semi- 
length L for a sample of Af identical loops follow from 

EM = n e n H (r + l)na 2 AfL 3 m n 2 £ 3 (62) 



where a is the loop aspect ratio (ratio between loop cross sectional diameter 
at the base and total length 2L) and r is the loop expansion factor. The 
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loop height for the important case of domina nt radiative losses follows to be 
llWhite et alJll986Uvan den Oord et alJll988l) 

^(i#) 1/S (^) 1/S( ^ r " 3(r + ir " 3 (63) 



A lower limit to H is found for r r ~ t c in the same treatment: 

^A^O^T 1 - (64) 



A/", a, and r are usually unknown and treated as free parameters within rea- 
sonable bounds. Generally, a small J\f is compatible with dominant radiative 
cooling. 



8.3 Quasi-static cooling loops 

Ivan den Oord fc Mewel l)l989h derived the energy equation of a cooling mag- 
netic loop in such a way that it is formally identical to a static loop 
<|Rosrier et al.lll 978h . by introducing a slowly varying flare heating rate that 
balances the total energy loss, and a possible constant heating rate during 
the flare decay. This specific solution thus proceeds through a sequence of 
different (quasi-)static loops with decreasing temperature. 

The general treatment involves continued heating that keeps the cooling 
loop at coronal temperatures. If this constant heating term is zero, one finds 
for free quasi-static cooling 

T(i) = r (l + </3r r ,o)~ 8/7 (65) 
L r (t) = L r , (l+V37v,or 4 (66) 
n e (t) = n ei0 (l + t/3r r! o)" 13/7 (67) 

where L r is the total radiative loss rate, and T r ^Q is the radiative loss time 
scale (I61II at the beginning of the flare decay. 

This prescription is equivalent to requiring a constant ratio between ra- 
diative and conductive loss times, i.e., in the approximation of T <; 20 MK 

(A oc T 1 / 4 ) 

T J113/4 

— = const w 0.18. (68) 

t c EM y ' 

Accordingly, the applicability of the quasi-static cooling approach can be 
supported or rejected based on the run of T and EM during the decay phase. 
Note, however, that a constant ratio (|68|l is not a sufficient condition to fully 
justify this approach. 



Stellar Coronae 41 



8.4 Cooling loops with continued heating 

Whether or not flaring loops indeed follow a quasi-static cooling path is best 
studied on a density-temperature diagram fFig. I18f) . Usually, characteristic 
values T = T a and n e = n e ^ as measured at the loop apex are used as 
diagnostics. For a magnetic loop in hydrostatic equilibrium, with constant 
cross section assumed, the RTV scaling law requires stable solutions 
(T,n e ) to be located where T 2 ps 7.6 x 10~ 7 n e L (for n e = ni). On a diagram 
of logT vs. logn e , all solutions are therefore located on a straight line with 
slope C = 0.5. FieurelTHlshows the path of a hydrodynamically simulated flare. 
The initial rapid heating leads to a rapid increase of T, inducing increased 
losses by conduction. As chromospheric evaporation grows, radiation helps 
to balance the heating energy input. The flare decay sets in once the heating 
rate drops. At this moment, depending on the amount of ongoing heating, 
the magnetic loop is too dense to be in equilibrium, and the radiative losses 
exceed the heating rate, resulting in a thermal instability. In the limit of no 
heating during the decay, that is, an abrupt turn-off of the heating at the 
flare peak, the slope of the path becomes 

din T T n „ 

C = = — =2 (69) 

ciln n e tt 

implying T(t) oc n^(t) = rv^(t) fsee lSerio et al.llT99ll for further discussion). 
Here, tt and r n are the e-folding decay times of the temperature and the elec- 
tron density, respectively, under the assumption of exponential decay laws. 
Only for a non- vanishing heating rate does the loop slowly recover and even- 
tually settle on a new equilibrium locus (Fig. IHfy . In contrast, if heating- 
continues and is very gradually reduced, the loop decays along the static so- 
lutions (^_^_0^] L _Obse^vationally, this path is often followed by large solar 
flares Uakimiec et alJ ll992). 



Applying this concept to stellar astronomy, we replace n e by the observ- 
able VEM and thus assume a constant flare volume, and further introduce 
the following generalization. In the freely cooling case after an abrupt heating 
turnoff, the entropy per particle at the loop apex decays on the thermody- 
namic decay time 

r td = 3.7x KT 4 -j^ [s] (70) 



wher e Tq is the flare temperature at the beginning of the decay ( Serio et al.l 



Il99l|) . When heating is p resent, we introduce a correction term F(() 
tlc /"Ttd l|Reale et al.lll997t tlc is the observed light curve decay time) 

r LC - 3.7 x lO- 4 -^^) [s] (71) 

F(() is therefore to be numerically calibrated for each X-ray telescope. With 
known F, l|71|l can be solved for L. This scheme thus offers i) an indirect 
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Fig. 18. Density-temperature diagram of a hydrodynamically simulated flare. 
The flare loop sta rts from an equilibrium (S-S, steady-state loop according to 
iRosner et aljfl973) : (a) and (b) refer to the heating phase; at (c), the heating 
is abruptly turned off, after which the loop cools rapidly (d, e), and only slowly re- 
covers tow^xda^nswequnibrium solution (f, g) due to constant background heating 
(from Ijakimiec et al.lFl992l) . 



method to study flaring loop geometries (L), ii) a way of determining the 
rate and decay time scale of continued heating via F (£) and T t( j, and iii) im- 
plications for the density decay time via r„ = £ t t- Conditions of applicability 
inclu de £ > 0.3 ( most values are found in the range of ~ 0.3 — l. lReale et al.l 
Il997|) and a resulting loop length L of less than one pressure scale height. 

8.5 Two-Ribbon flare models 

An approach that is entirely based on continuous heating (as opposed to cool- 
ing) was developed for the two-ribbon (2-R) class of solar flares. An example 
of this flare type is shown in Fig. ^| The 2-R flare model devised initially 
bv lKopp &: Polettcl l|l984jl is a parameterized magnetic-energy release model. 
The time development of the flare light-curve is completely determined by the 
amount of energy available in non-potential magnetic fields, and by the rate 
of energy release as a function of time and geometry as the fields reconnect 
and relax to a potential-field configuration. It is assumed that a portion of the 
total energy is radiated into the observed X-ray band, while the remaining 
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energy will be lost by other mechanisms. 2-R flares are well established for 
the Sun fFig. 119( 1: they often lead to large, long-duration flares that may be 
accompanied by mass ejections. 




Fig. 19. Trace image of a flaring magnetic loop arcade. 



The magnetic fields are, for convenience, described along meridional 
planes on the star by Legendre polynomials P n of order n, up to the height 
of the neutral point; above this level, the field is directed radially, that is, the 
field lines are "open" . As time proceeds, field lines nearest to the neutral line 
move inward at coronal levels and reconnect at progressively larger heights 
above the neutral line. The reconnection point thus moves upward as the 
flare proceeds, leaving closed magnetic-loop systems underneath. One loop 
arcade thus corresponds to one N-S aligned lobe between two zeros of P n in 
latitude, axisymmetrically continued over some longitude in E-W direction. 
The propagation of the neutral point in height, y(t), with a time constant to, 
is prescribed by (y in units of i?*, measured from the star's center) 

y(t) = 1 + |= (l - e"*/*) (72) 

H(t) = [y(t) - 1}R* (73) 

and the total energy release of the reconnecting arcade per radian in longitude 
is equal to the magnetic energy lost by reconnection, 
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2±. = J_2n(n + l)(2n + l) 2 R 3 .B 2 I 12 (n) [ V - } y . . ^ (74) 
dy 8?r v A ; * v ; [n + (n + l)y 2n+1 ] 3 v ; 



(75) 



dE _ dE dy 
dt dy dt 

()Poletto et alJll988^ . In G2J), H m is the maximum height of the neutral point 
for t —>■ oo; typically, H m is assumed to be equal to the latitudinal extent of 
the loops, i.e., 

H m ps -^-j- R* (76) 
n + 1/2 

for n > 2 and H m — (tv/2)R* for n = 2. Here, i? is the surface magnetic 
field strength at the axis of symmetry, and i?* is the stellar radius. Finally, 
Ii2(n) corresponds to J[P n {cos9)] 2 d(cos9) evaluated between the latitudinal 
borders of the lobe (zeros of dP n /d#), and 9 is the co-latitude. 

The free parameters are B and the efficiency of the energy-to-radiation 
conversion, q, both of which determine the normalization of the light curve; 
the time scale of the reconnection process, to, and the polynomial degree n 
determine the duration of the flare; and the geometry of the flare is fixed 
by n and therefore the asymptotic height H m of the reconnection point. The 
largest realistic 2-R flare model is based on the Legendre polynomial of degree 
n = 2; the loop arcade then stretches out between the equator and the stellar 
poles. Usually, solutions can be found for many larger n as well. However, 
because a larger n requires larger surface magnetic field strengths, a natural 
limit is set to n within the framework of this model. Once the model solution 
has been established, further parameters, in particular the electron density 
n e , can be inferred. 



8.6 A magnetohydrodynamic model 

Some scaling laws have been obtaine d from simulations based on the fu ll set 
of magnetohydrodynamic equations l|flhiha,ta fc Yokovamal Fmflfll 120021) . For 
the flare peak temperature T, the loop magnetic field strength B, the pre- 
flare loop electron density no, and the loop semi-length L, one finds, under 
the condition of dominant conductive cooling (appropriate for the early phase 
of a flare), 

T w 1.8 x 10 4 B 6 / 7 < 1,7 L 2 ' 7 [K]. (77) 

The law follows from the balance between conduction cooling (ex T 7 / 2 /£ 2 . I58|) 
and magnetic reconnection heating (oc B 3 /L). Assuming loop filling through 
chromospheric evaporation and balance between thermal and magnetic pres- 
sure in the loop, two further "pressure-balance scaling laws" follow: 

EM fa 3 x lO- 17 B- 5 n 3 /2 T 17/2 [cm" 3 ] (78) 
EM ps 2 x 10 8 L 5/3 no /3 T 8/3 [cm" 3 ]. (79) 
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An alternative scaling law applies if the density development in the initial 
flare phase is assumed to follow balance between evaporation enthalpy-flux 
and conduction flux, although the observational support is weaker, 

EM « f x f0" 5 B" 3 7io /2 T 15/2 [crrT 3 ]. (80) 

And third, a steady solution is found for which the radiative losses balance 
conductive losses. This scaling law applies to a steady loop, 

Jf0 13 T 4 i [cm- 3 ] for T < 1 7 K 
tM ~ \ 1 20 T 3 L [cm- 3 ] for T > 1 7 K i8iJ 

and is equivalent to the RTV scaling law l|33fl . 

The advantage of these scaling laws is that they make use exclusively of 
the flare-peak parameters T, EM, B (and the pre-flare density n ) and do 
not require knowledge of the time evolution of these parameters. 



8.7 Observations of stellar X-ray flares 

One of the main results that have come from extensive modeling of stellar 
X-ray flares is that extremely large stellar flares require large volumes under 
all realistic assumptions for the flare density, i.e., flaring complexes of mag- 
netic loops must either be high, or must spread across a large surface area. 
This is because, first, the energy derives from the non-potential portion of 
the magnetic fields that are probably no stronger than a few fOO G in the 
corona; and second, small-loop models require higher pressure to produce the 
observed luminosity, hence requiring excessively strong magnetic fields. 

There is extensive literature discussing individual flares observed on a va- 
riety of stars. 1 will not discuss these results indi vidually here. The interested 
reader may consult the compilation of results in lGiideJ l|2004j) and the refer- 
ences given therein. This section summarizes some systematic trends found 
in stellar flares. 

When the flare energy release evaporates plasma into the corona, heating 
and cooling effects compete simultaneously, depending on the density and 
temperature profiles in a given flare, ft is therefore quite surprising to find a 
broad correlation between peak temperature T p and peak emission meas ure 
EM p , as illustrated in Fig. HOI for the sample reported in iGiidell l)2004|) . A 
regression fit gives (for 66 entries) 

EM p cx Tp 4 - 30±0 - 35 . (82) 

The correlation overall indicates that larger flares are hotter. A s imilar rela- 
tion was reported previously for solar flares ^Feldman et aljfl995l) . 

It is interesting to note that this correlation is similar to the T — Lx 
correlation for the "non-flaring" coronal stars in Fig.0at cooler temperatures. 
This same sample is plotted as triangles in Fig. 1201 again only for the hotter 
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Fig. 20. Le/t: Peak temperatures and EMs of the flares ('from lOudel I2004I and 
references therein). Key to the symbols: Filled circles: XMM-Newton observations. 
Open circles: ASCA or BeppoSAX observations. Small diamonds: observations from 
other satellites. The solid line shows a regression fit 1821 . Triangles represent non- 
flaring parameters of G stars, referring to the hotter plasma component in 2-T 
spectral fits to ROSAT data. - Right: Theoretical EM-T relations based on the 
reconnection model by Shibata & Yokoyama, showing lines of constant loop length 
L and lines of constant magnetic field strength B. Hatched areas are loci reported 
for solar flares, and other symbols refer to individual stellar flares in star-forming 
regio ns (figure courtesy of K. Shibata and T. Yokoyama, after lShibata fc Yokovamal 



plasma component. The stars follow approximately the same slope as the 
flares, albeit at cooler temperatures, and for a given temperature, the EM is 
higher. This trend may suggest that flares systematically contribute to the 
hot plasma component, although we have not temporally averaged the flare 
temperature and EM for this simple comparison. 

In the context of the magnetohydrodynamic scaling laws presented in 
1771) (|78|l , the observed loci of the flares require loop magnetic field strengths 
similar to solar flare values (B w 10 — 150 G) but the loop lengths must 
increase toward larger flares. This is seen in Fig.[2U]where lines of constant L 
and B are plotted for this flare model. Typical loop lengths are thus of order 
L m 10 11 cm in this interpretation. 

8.8 The "Neupert Effect" 

In the framework of chromospheric evaporation, we assume that the emitted 
radio emission or its flux at Earth, Pr(4) is at any time proportional to the 
deposition rate of kinetic energy by nonthermal electrons into the plasma of 
the chromosphere. This plasma is thereby heated and "evaporates" into the 
corona. We define the conversion factor a as the ratio between the thermal 
energy flux being deposited in the chromosphere by nonthermal electrons 
and the observed radio flux density. Conversion into other forms of energy 
(mechanical, turbulence, etc.) occurs in parallel. Our simplification consists, 
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Fig. 21. Flare on Proxima Centauri, observed with XMM-Newton. The top panel 
shows the X-ray light curve and the much shorter U band flare (around 17 UT). 
The bottom panel shows the Ovu He-like triplets observed during various time 
intervals of the flare. The locations of the r, i, and / lines are marked by vertical 
lines. The resulting electron densities are given in the top panel by the crosses, 
where the horizontal arm lengths indicate the time intervals over wh ich the data 
were i ntegrated, and the right axis gives the logarithmic scale f after iGiidel et all 
20023). 



first, in the assumption that these other energy transformations do not inter- 
act so that a remains constant over the course of the relatively short radio 
flare; second, we will assume that the dominant losses of the thermal energy 
occur via radiation; we will neglect energy loss into cool, non- X-ray emitting 
plasma, e.g. by conduction, and cooling by adiabatic expansion. Also, we 
assume that there is no direct heating in parallel to the chromospheric evap- 
oration. Third, we will use one temperature parameter for the flare plasma 
at any given time: The temperature T can be considered as describing an 
isothermal plasma dominating the total losses from the flaring loops. 

The rate of change of the total thermal energy E in a plasma of volume 
V with electron density n e is determined by the influx of kinetic energy and 
by radiation, hence the energy conservation equation for the thermal plasma 
is 
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j t {Zn e kTV) = aF R (t) - n\V A{T) (83) 

where A(T) is the total luminosity of a plasma with unit emission measure 
(EM) at a temperature of T. Note again that 1)83(1 does, by definition of a, not 
describe the total energy budget of the flare, but merely the energy conversion 
of interest here. For large flare temperatures (3> 20 MK), A is dominated by 
bremsstrahlung losses, roughly scaling as T 1 / 2 . For somewhat more moderate 
temperatures ( > 10—20 MK), losses via line emission become important, 
with the (isothermal) approximation for A 

A(T) = 1.86 • l(T 25 T 1/4 [erg s" 1 cm 3 ] (84) 

(see after equation 1591- for simplicity, we define here EM = n\V with nn = 
0.85n e and a corresponding correction in A). In the absence of kinetic energy 
influx aF R (t), the thermal energy decays radiatively (neglecting conduction), 
and hence we define an e-folding decay time r for the thermal energy as a 
function of the two independent variables temperature T and density n e at 
a given instant, 



dE{t) 



dt 

so that 



F R =0 



T 



E(t) (85) 



(n e (t),T(t)) 



. . „, NN E(t) 3kT 
r(n„(t), T(t)) = — = 86 

where L ra d is the total luminosity (approximately = X-ray luminosity, Lx, 
but also contributions from lower temperatures, e.g., in the ultraviolet, which 
we neglect here). Then, 

dECt) , , , , , , E(t) , s 

= aF R (t) - L rad (i) = aF R (t) - (87) 

The general solution of the inhomogeneous, linear differential equation l|87|l 
reads 

E(t) = e Jt o I Eq + a t R {u)e Jt o du I 

where Eq = E(to) for a fixed to before the flare start. The integration over r _1 
is along the time axis. If the initial thermal energy content can be neglected 
(E Q = 0), then 

E(t) = a f F R (u)e- £ T[vrldy du. 



We define an average decay constant for any time interval [u, t] by 



r 



-l 



\u,t) = S » T{y) ~ ldy . (90) 
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Then 

E(t) = a [ F R {u)e- {t - u)/f ^ u) du, (91) 

Jto 

i.e., for constant r the energy profile is the convolution of the kinetic energy 
influx with an exponential function. Equation H91(l is a generalized form of 
the Neupert effect. In the limit t — > oo, (|87(l and (|fJ 1|) become 

dE(t) 

= afk(«) (92) 

= ot f F R (u)du, (93) 

i.e., the total energy content of the plasma is the integral of the kinetic energy 
influx (and radiation is inhibited). Equations 192(1 and 1931) remind us of the 
classical formulation of the Neupert effect, with the observed X-ray losses 
replaced by the thermal energy content of the plasma. These two equations 
are applicable only for the increasing portion of the soft X-ray light curve 
or, correspondingly, for the time interval where F R ^ 0. On using l|86[) . we 
obtain the generalized Neupert effect for the light curve, 

L rad (i) = " / F R (u)e-^/ f ^du. (94) 

Note the importance of the thermal energy decay time t = 3kT / (n e A(T)) . 
Serio et al. (1991) and Jakimiec et al. (1992) find Tn~ 2 = const for the radia- 
tive cooling phase. If the temperatures are very high and the bremsstrahlung 
approximation applies, then r = T 1//2 rt -1 and thus r = const in (|94|l . How- 
ever, in the case of more moderate temperatures (I84|) . only for T 3 / 4 n _1 = 
const does the classical Neupert effect approximately apply to the light curve. 
The actual functional dependence ofip(T) is more complicated. Thus, in gen- 
eral, the relevant parameters to be investigated for the Neupert effect are the 
energies. Further, I/ ra d(i) describes all radiative losses across the electromag- 
netic spectrum; the observed X-ray luminosity Lx generally depends on the 
selected bandpass and therefore constitutes only a lower limit to £ ra d- 

Despite these caveats, it is surprising that the Neupert effect is often well 
seen, at least qualitatively, in solar radio and X-ray light curves. Instead of 
radio emission, bursts can be monitored in the U band: because they are likely 
to be a prompt reaction to the bombardment of the chromosphere by the same 
electron population that indices gyrosynchrotron emission, the Neupert effect 
should hold as well. A stellar example is shown in Fig. |^ These observations 
give clear evidence that at least some giant flares on active stars are subject 
to similar evaporation physics as known from the solar corona. 
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9 The statistics of flares 

The study of coronal structure confronts us with several problems that are 
difficult to explain by scaling of solar coronal structure: i) Characteristic coro- 
nal temperatures increase with increasing magnetic activity, ii) Characteristic 
coronal densities are typically higher in active than in inactive stars, and pres- 
sures in hot loops can be exceedingly high, iii) The maximum stellar X-ray 
luminosities exceed the levels expected from complete coverage of the surface 
with solar-like active regions by up to an order of magnitude, iv) Radio obser- 
vations reveal a persistent population of non-thermal high-energy electrons 
in magnetically active stars even if the lifetime of such a population should 
only be tens of minutes to about an hour under ideal trapping conditions in 
coronal loops and perhaps much less due to efficient scattering of electrons 
into the chromosphere. Several of these features are reminiscent of flaring, 
as are some structural elements in stellar coronae. If flares are important for 
any of the above stellar coronal properties indeed, then we must consider the 
effects of frequent flares that may be unresolved in our observations but that 
may make up part, if not all, of the "quiescent" emission. 

9.1 Stochastic variability - what is "quiescent emission"? 

The problem has been attacked in several dedicated statistical studies. Early 
statistical investigations (fluctuation analysis in light curves) remained am- 
biguous, reporting a significant amount of variability, although not necessarily 
being due to flar es, or statistical absence of low-level flaring within the sensi- 
tivity limits fe.g-. lAmbruster et alJl987tlCollura et alJl988ilPallavicini et all 
1990). There are indications in newer data that M dwarfs are continuously 
variable on short ( 1 day) time scales, and that the luminosity distribution 
is very similar to the equivalent distribution derived for solar flares, which 
suggests that the overall stellar light curves of dM stars are variab le in the 
same way as a statistical sample of solar flares l|Marino et al.ll200ol) . 

Very long light curves obtained from the EUVE satellite reveal an as- 
tonishing level of continuous variability in active M dwarfs. Some of those 
data can be used to investigate statistical properties of the occurrence rate 
of flares as a function of total emitted energy. The increased sensitivity of 
XMM-Newton and Chandra is now revealing extreme levels of activity. Some 
X-ray light curves show no steady time interval exceeding a few tens of min- 
utes within the sensitivity limit. In the day-long light curve in Fig. 1231 no 
more than 30%, and probably much less, of the average X-ray emission of 
UV Cet can be attributed to any sort of steady emission, even outside the 
obvious, large flares. On the contrary, almost the entire light curve is resolved 
into frequent, stochastically occurring flares of various amplitudes. 

In the solar corona, the flare rate increases steeply towar d lower radiative 
energ ies, with no evidence (yet) for a lower threshold (e.g., iKrucker fc Bend 
1998). Figure El shows an example of a GOES light curve in the 1.5-12 keV 
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Fig. 22. A long light curve of the dMe star AD Leo, obtained by the D S instrument 
on EU VE. Most of the discernible variability is due to flares (after IGiidel et alJ 
l2003al) . 



Chandra LETGS 0th order 
UV Celi B 




Fig. 23. Light curve of UV Ceti B, observed with the Chandra LETGS/HRC 
durin g about 1 day. Not e the logarithmic flux axis (figure courtesy of M. Audard, 
after lAudard et alJ l2003 ) . 



range, purposely selected during an extremely active period in November 
2003. While the GOES band is harder than typical bands used for stellar 
observations, it more clearly reveals the level of the underlying variability (a 
typical detector used for stellar observations would see much less contrast). 
If the solar analogy has any merit in interpreting stellar coronal X-rays, then 
low- level emission in stars that do show flares cannot be truly quiescent, that 
is, constant or slowly varying exclusively due to long-term evolution of active 
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Fig. 24. GOES full-disk solar X-ray light curve, observed in the 1.5-12 keV band 
in November 2003. The abscissa gives time after 2003 November 1, 7:12 UT in days. 

regions, or due to rotational modulation. A measure of flare rates is therefore 
not meaningful unless it refers to flares above a given luminosity or energy 
threshold. This is - emphatically - not to say that steady emission is absent 
in magnetically active stars. However, once we accept the solar analogy as a 
working principle, the question is not so much about the presence of large 
numbers of flares, but to what extent they contribute to the overall X-ray 
emission from coronae. 



9.2 The flare-energy distribution 

The suggestion that stochastically occurring flares may be largely responsible 
for coronal heati ng is known as the "microflare" or "nanoflare" hypothesis 
in solar physics l|Parkerlll988lh Observationally, it is supported by evidence 
for the presence of nume rous small-scale flare events occurring in the solar 
corona at any time fe.g.. lLin et al.lfl984() . Their distribution in energy is a 
power law, 

clN 

% = < 95 > 

where dN is the number of flares per unit time with a total energy in the 
interval [E, E + dE], and k is a constant. If a > 2, then the energy integration 
(for a given time interval) diverges for i? m i n — > 0, that is, by extrapolating the 
power law to sufficiently small flare energies, any energy release power can be 
attained. This is not the case for a < 2. Solar st udies have repeated ly resulted 
in a values of 1.6 — 1.8 for ordinary solar flares llCrosbv et alj|l993 |). but some 
recent studies of low-l evel flaring suggest a = 2.0—2.6 l|Krucker fc BenJl998t 
IParnell fc Jiir^OOOh . 

Relevant stellar studies have been rare (see Table 0J. Early investigations 
lumped several stars together to produce meaningful statistics. Full forward 
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modeling of a supe rposition of stochastic flar es was applied to EU V and X- 



rav light curves bvlKashvaD et al. 




2002) andlGiidel et all J2003al) based on 


Monte Carlo simulations, and by 


Arzner & Gudell 


2004 


based on an ana- 



lytical formulation. The results of these investigations are in full agreement, 
converging to a w 2.0 — 2.5 for M dwarfs (Table 0). If the power-law flare 
energy distribution extends by about 1-2 orders of magnitude below the ac- 
tual detection limit in the light curves, then the entire emission could be 
explained by stochastic flares. The coronal heating process in magnetically 
active stars would - in this extreme limit - be one solely due to time- dependent 
heating by flares, or, in other words, the X-ray corona would be an entirely 
hydrodynamic phenomenon rather than an ensemble of hydrostatic loops. 




Fig. 25. Left: The rate of flares above a threshold of 10 erg in total radiated 
X-ray energy is plotted against the low-level luminosity for several stars, together 
with a regression fit. Right: Flare energy distribution for AD Leo, using a flare 
identification algori thm for an obs ervation with EUVE (both figures courtesy of M. 
Audard, after lAudard et al.ll200ol) . 



Table 3. Stellar radiative flare-energy distributions 



Star sample 


Photon energies 


log (Flare 


Q 


References 




[keV] 


ener 


?ies) a 






M dwarfs 


0.05-2 


30.6 


-33.2 


1.52±0.08 Collura et al. (1988) 


M dwarfs 


0.05-2 


30.5 


- 34.0 


1.7L0.1 


Pallavicini et al. (1990) 


RS CVn binaries EUV 


32.9 


- 34.6 


1.6 


Osten fc Brown (1999) 


Two G dwarfs 


EUV 


33.5 


- 34.8 


2.0-2.2 


Audard et al. 1999) 


F-M dwarfs 


EUV 


30.6 


- 35.0 


1.8-2.3 


Audard et al. (2000) 


Three M dwarfs 


EUV 


29.0 


- 33.7 


2.2-2.7 


KashvaD et al. C2002) 


AD Leo 


EUV&0.1-10 


31.1 


- 33.7 


2.0-2.5 


Giidel et al. f 2003a) 


AD Leo 


EUV 


31.1 


- 33.7 


2.3 ± 0.1 


Arzner fc Giidel (2004) 



"Total flare-radiated X-ray energies used for the analysis (in ergs). 
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Fig. 26. Light curves of HR 1099 (left) and UX Ari (right) obtained in the two 
senses of circular polarization. The brighter of the two polarized fluxes varies rapidly 
and has been interpreted as 100% polarized coherent emission superimposed on a 
gradually changing gyrosynchrotron component (White & Franciosini 1995; figures 
from S. M. White). 

9.3 Microflaring at radio wavelengths 

Quiescent radio emission can apparently persist for quite long periods. Losses 
by collisions (|llfl require a very low ambient electron density to maintain 
the electron population. Alternatively, electrons could be frequently injected 
at ma ny coronal sites. Based on spectral observations. IWhite fc Franciosinil 
(1995) suggest that the emission around 1.4 GHz shown in Fig.lSrjIis composed 
of a steady, weakly polarized broad-band gyrosynchrotron component plus 
superimposed, strongly and oppositely polarized, fluctuating plasma emission 
that is perceived as quasi-steady but that may occasionally evolve into strong, 
polarized flare emission. Continual flaring may thus also reflect in radio light 
curves. 

10 A flare-heating approach 

The quiescent emission of synchrotron radiation requires the continuous pres- 
ence of nonthermal electrons, while X-rays are emitted by the thermal bulk 
of typically 10 6 to a few times 10 7 K. For optically thin radio emission (v > 5 
GHz), the radio luminosity is 

L R = Att^Vr [erg s" 1 Hz" 1 ] , (96) 



Stellar Coronae 55 



where Vr is the radio source volume. The density of nonthermal electrons 
will be assumed to be distributed in energy according to a power law, 

n(e) = ^— ^ N (—) [cm" 3 erg" 1 ]. (97) 

lower cutoff at £q w 10 keV = 1.6 • 10 -8 erg is compatible with most accelera- 
tion processes proposed for stellar coronae and with solar flare observations. 
We keep £o fixed at 10 keV in the following considerations. It is not a sensi- 
tive value for synchrotron emission since the latter becomes appreciable only 
at higher energies. We assume a homogeneous sourc e for s i mplici ty. The gy- 
rosynchrotron emissivity is approximately given by iDulkl l|l985l) (sum of x 
and o mode, see equation I46fl 

f \ 1.22-0.905 

r, « 1.8- 3.3 x 10- 24 BN lO-°- 525 (sin0)-°- 43+0 - 655 ( J (98) 

= 1.8 • d{B, 5, v, 6) N [erg s^Hz" 1 cm" 3 sterad" 1 ] , (99) 

where vb is the electron gyrofrequency and 8 is the angle between the line 
of sight and the magnetic field. In a steady-state situation, the density of 
nonthermal electrons is given by 

n{e) = ^r(e) (100) 

where n; n (e) is the total number of electrons of energy e accelerated per unit 
time, and t(e) is the electron lifetime. Let 

n in ( e ) = ho,in C^j (101) 



r(e)=T 1-1 (102) 

With (|100l) . the power law index of the electrons (|97|) is S i=s k — a. 

Let a be the fraction of the energy that goes into accelerated particles, 
and b the fraction of the total coronal energy ultimately radiated into the 
observed X-ray band. Since some of the thermal energy is lost by conduction 
and other processes, b < 1. Then, X-rays are related to the total energy input 

1 f 00 1 
E=- h in (e)e de = -L X - (103) 

a J SO " 

Using (|9*rj|) - (| 103(1 . the relation between Lr and Lx becomes 

a _i / a — 1 



L R = 1.8 • 47rt?(fl, u, 9, 5) - £o L t I — - + 1 ) L x (104) 
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where we require 6 > 1 and a > 2 — 5 for convergence. Equation I|1U4|) is 
general and includes different possible scenarios. Let us select typical param- 
eters for stellar observations, viz. 6 = 3 (implying a > — f), v = 5 GHz, and 
6 = 30°; then 

L R = 3.5 • f(T 22 5 2 - 48 | r (a + f ) L x (105) 

Late- type main-sequence and subgiant stars appear to follow a linear relation 
between Lr and Lx- L \ ~ 10 15,5 I/r (the co efficient is 0.5-1 dex smaller for 
RS CVn-type binaries: iGiidel fe Benzlll993l see Fig. 127)1 . We thus suggest 
numerical relations between a, &, r, and B for active stars: 

R2.48 a f 9.0 • 10 s dMe, dKe, BY Dra , , 

6 T0la + j ~ 1 5.4 • 10 6 RS CVn, Algols, PTTS, FK Com. [ Ub) 



Two scenarios are possible: If the acceleration efficiency a is close to unity, 
the nonthermal electrons first emit a small fraction of the total kinetic energy 
as synchrotron radiation before they lose most of their energy by collisions, 
thereby heating the X-ray emitting corona ("causal relation" between non- 
thermal and thermal energy). It is unknown whether fully efficient accelera- 
tors are realized in nature. 

The other scenario implies a common energy release, most of it heating 
the corona by thermal processes (e.g. Ohmic heating), and a much smaller 
fraction a accelerating particles. Some of the latter energy may also end up 
as heat, but this is negligible. Therefore, the thermal plasma and nonthermal 
electrons radiate independently ("common origin" scenario). This scenario 
can explain the relation found by assuming that coronal heating in active 
stars necessarily implies electron acceleration to relativistic energies. The 
fraction between thermal and nonthermal energy release is the crucial, but 
unknown number. It is unclear whether it can be constant over many different 
types of stars. 

The lifetime r (or To and a) which is relevant here for the gyrosynchrotron 
emitting high-energy particles, depends on the process that scatters ener- 
getic electrons into the loss-cone of the trapped particle velocity distribution, 
such as collisions, whistler wave instability, or cyclotron maser action. Their 
effects on the relation between trapping time and particle energy are op- 
posite. As an example, let us assume r(e) = const (i.e. a = 0) and use 
solar active region and flare values for B (w 100 G) and a/b (order of unity; 
IPennisI Il988t IWithbroe fc Noveslll977i Equation (fT77Hf> then would require 
r « 10 s for dM(e), dK(e), and BY Dra stars. The lifetime is an average 
over the time of flight of the particles that are lost immediately (moving 
parallel to the field lines) and the trapped particles temporarily residing in 
a lower density plasma. The estimated values for r are compatible with the 
minimum observed time scale of variations in the "quiescent" emission (e.g. 
I.Tackson et al J 1 9891. F or the halo of RS CVn binaries, B can be as low as w 10 
G llMutel et al.lll985h . Then, from dUl, t is of order of many hours, com - 
patible with observed long time scales Ce.g.. lMassi fc Chiuderi-Dragd ll992). 
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Fig. 27. C orrelation be tween radio and X-ray luminosities for magnetically active 
stars (from lGiidelll2002h . 
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